MMR-and-RMR-workflow-Clutch59
Elliott Schmidt
2024-06-26
Load libraries (packages)
System1 - Dell
Importing data from firesting for resting
preexperiment_date <- "04 April 2023 12 20PM/All"
postexperiment_date <- "04 April 2023 06 02PM/All"
##--- last fish run in trial ---##
experiment_date <- "04 April 2023 02 58PM/Oxygen"
experiment_date2 <- "04 April 2023 02 58PM/All"
firesting <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_date,"data raw/Firesting.txt"),
delim = "\t", escape_double = FALSE,
col_types = cols(`Time (HH:MM:SS)` = col_time(format = "%H:%M:%S"),
`Time (s)` = col_number(), Ch1...5 = col_number(),
Ch2...6 = col_number(), Ch3...7 = col_number(),
Ch4...8 = col_number()), trim_ws = TRUE,
skip = 19)
Cycle_1 <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_date2,"slopes/Cycle_1.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE)
Cycle_last <-read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_date2,"slopes/Cycle_20.txt"), # custom
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE) System2 - Asus
Importing data from firesting for resting
preexperiment_date_asus <- "04 April 2023 12 03PM/All"
postexperiment_date_asus <- "04 April 2023 05 31PM/All"
##--- last fish run in trial ---##
experiment_date_asus <- "04 April 2023 02 10PM/Oxygen"
experiment_date2_asus <- "04 April 2023 02 10PM/All"
firesting_asus <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_date_asus,"data raw/Firesting.txt"),
delim = "\t", escape_double = FALSE,
col_types = cols(`Time (HH:MM:SS)` = col_time(format = "%H:%M:%S"),
`Time (s)` = col_number(), Ch1...5 = col_number(),
Ch2...6 = col_number(), Ch3...7 = col_number(),
Ch4...8 = col_number()), trim_ws = TRUE,
skip = 19)
Cycle_1_asus <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_date2_asus,"slopes/Cycle_1.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE)
Cycle_last_asus <-read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_date2_asus,"slopes/Cycle_20.txt"), # custom
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE) Chamber volumes
chamber1_dell = 0.04650
chamber2_dell = 0.04593
chamber3_dell = 0.04977
chamber4_dell = 0.04860
chamber1_asus = 0.04565
chamber2_asus = 0.04573+0.00385
chamber3_asus = 0.04551+0.00322
chamber4_asus = 0.04791+0.00277
Date_tested="2023-04-04"
Clutch = "59"
Male = "CARL354"
Female = "CARL355"
Population = "Arlington reef"
Tank =281
salinity =36
Date_analysed = Sys.Date() Replicates
1
Enter specimen data
Replicate = 1
mass = 0.0006711
chamber = "ch4"
Swim = "good/good"
chamber_vol = chamber4_dell
system1 = "Dell"
Notes="check max"
##--- time of trail ---##
experiment_mmr_date <- "04 April 2023 02 26PM/Oxygen"
experiment_mmr_date2 <- "04 April 2023 02 26PM/All"
firesting_mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_mmr_date,"data raw/Firesting.txt"),
delim = "\t", escape_double = FALSE,
col_types = cols(`Time (HH:MM:SS)` = col_time(format = "%H:%M:%S"),
`Time (s)` = col_number(), Ch1...5 = col_number(),
Ch2...6 = col_number(), Ch3...7 = col_number(),
Ch4...8 = col_number()), trim_ws = TRUE,
skip = 19) ## New names:
## • `Ch1` -> `Ch1...5`
## • `Ch2` -> `Ch2...6`
## • `Ch3` -> `Ch3...7`
## • `Ch4` -> `Ch4...8`
## • `Ch 1` -> `Ch 1...9`
## • `Ch 2` -> `Ch 2...10`
## • `Ch 3` -> `Ch 3...11`
## • `Ch 4` -> `Ch 4...12`
## • `('C)` -> `('C)...15`
## • `('C)` -> `('C)...16`
## • `Ch 1` -> `Ch 1...18`
## • `Ch 2` -> `Ch 2...19`
## • `Ch 3` -> `Ch 3...20`
## • `Ch 4` -> `Ch 4...21`
## • `Ch1` -> `Ch1...22`
## • `Ch2` -> `Ch2...23`
## • `Ch3` -> `Ch3...24`
## • `Ch4` -> `Ch4...25`
## • `Ch1` -> `Ch1...26`
## • `Ch2` -> `Ch2...27`
## • `Ch3` -> `Ch3...28`
## • `Ch4` -> `Ch4...29`
## • `` -> `...31`## Warning: One or more parsing issues, call `problems()` on your data frame for details,
## e.g.:
## dat <- vroom(...)
## problems(dat)Cycle_1.mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_mmr_date2,"slopes/Cycle_1.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE) ## New names:
## • `` -> `...8`Background rates
Pre-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",preexperiment_date,"slopes"))
pre_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_pre1 <- pre_cycle1 %>% calc_rate.bg()##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] 0.0004003393post-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",postexperiment_date,"slopes"))
post_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_post1 <- post_cycle1 %>% calc_rate.bg() ##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] 0.001184524Resting metabolic rate
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch4
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 0.94 3.86
## -----------------------------------------
#### subset data
Tstart.row=which(firesting2$TIME ==Cycle_1$Time[1], firesting$TIME)
Tstart.dTIME=as.numeric(firesting2[Tstart.row[1], "dTIME"]) # custom
Tend.row=which(firesting2$TIME ==tail(Cycle_last$Time, n=1), firesting$TIME)
Tend.dTIME=as.numeric(firesting2[Tend.row, "dTIME"])
apoly_insp <- firesting2 |>
subset_data(from=Tstart.dTIME,
to=Tend.dTIME,
by="time")
apoly_insp <- inspect(apoly_insp, time=1, oxygen=2)##
## # print.inspect # -----------------------
## dTIME ch4
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 7 8 9 10 11 12 14 16 18 19 20 21 22 23 24
## Minimum and Maximum intervals in uneven Time data:
## [1] 0.94 3.86
## -----------------------------------------
Extract rates
apoly_cr.int <- calc_rate.int(apoly_insp,
starts=(195+45+300),
wait=45,
measure=255,
by="time",
plot=TRUE) ##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from all replicates ...
## -----------------------------------------##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from selected replicates...
## To plot others modify 'pos' input.
## -----------------------------------------adjust rates for background
apoly_cr.int_adj <- adjust_rate(apoly_cr.int,
by = bg_pre,
by2 = bg_post,
time_by = Tstart.row[1], # custom
time_by2 = Tend.row,
method = "linear")## Warning: adjust_rate: One or more of the timestamps for the rate(s) in 'x' do not lie between the timestamps for the 'by' and 'by2' background rates.
## Ensure this is correct. The adjustment value has been calculated regardless by extrapolating outside the background rates time window.## adjust_rate: Rate adjustments applied using "linear" method.Converting units
apoly_cr.int_adj2 <- apoly_cr.int_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253) ## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.Plot curve
ggplot(as.data.frame(apoly_cr.int_adj2$summary), aes(x=row, y=rate.output*-1)) +
geom_point() +
stat_smooth(method = "lm", formula = y~poly(x, 2), color="red") +
theme_classic()
Rate filtering
apoly_rmr <- apoly_cr.int_adj2 |>
select_rate(method ="rsq", n=c(0.95,1)) |>
select_rate(method="lowest", n=6) |>
plot(type="full") |>
summary(export = TRUE)## select_rate: Selecting rates with rsq values between 0.95 and 1...## ----- Selection complete. 0 rate(s) removed, 20 rate(s) remaining -----## select_rate: Selecting lowest 6 *absolute* rate values...## ----- Selection complete. 14 rate(s) removed, 6 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: 10 1 236.5052 -0.009597093 0.979 NA 4840 5095 14416.69
## 2: 11 1 244.9909 -0.009817408 0.965 NA 5380 5635 14956.72
## 3: 15 1 253.3372 -0.009070945 0.969 NA 7540 7795 17116.78
## 4: 16 1 263.8381 -0.009387837 0.960 NA 8080 8335 17656.71
## 5: 18 1 271.2200 -0.009233169 0.970 NA 9160 9415 18736.69
## 6: 20 1 283.9111 -0.009366765 0.955 NA 10240 10495 19816.69
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 14671.69 98.088 95.597 -0.009597093 0.0008005135 -0.01039761 -0.01039761
## 2: 15211.69 97.856 95.593 -0.009817408 0.0008408518 -0.01065826 -0.01065826
## 3: 17371.72 97.949 95.551 -0.009070945 0.0010022040 -0.01007315 -0.01007315
## 4: 17911.79 97.792 95.674 -0.009387837 0.0010425412 -0.01043038 -0.01043038
## 5: 18991.70 97.988 95.616 -0.009233169 0.0011232115 -0.01035638 -0.01035638
## 6: 20071.70 97.933 95.604 -0.009366765 0.0012038859 -0.01057065 -0.01057065
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.0486 0.0006711 NA 36 27 1.013253 -0.1180957
## 2: %Air sec 0.0486 0.0006711 NA 36 27 1.013253 -0.1210562
## 3: %Air sec 0.0486 0.0006711 NA 36 27 1.013253 -0.1144105
## 4: %Air sec 0.0486 0.0006711 NA 36 27 1.013253 -0.1184679
## 5: %Air sec 0.0486 0.0006711 NA 36 27 1.013253 -0.1176275
## 6: %Air sec 0.0486 0.0006711 NA 36 27 1.013253 -0.1200611
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -175.9733 NA mgO2/hr/kg -175.9733
## 2: -180.3847 NA mgO2/hr/kg -180.3847
## 3: -170.4821 NA mgO2/hr/kg -170.4821
## 4: -176.5280 NA mgO2/hr/kg -176.5280
## 5: -175.2756 NA mgO2/hr/kg -175.2756
## 6: -178.9020 NA mgO2/hr/kg -178.9020
## -----------------------------------------Results
results <- data.frame(Clutch = Clutch,
Replicate =Replicate,
Male=Male,
Female=Female,
Population = Population,
Tank = Tank,
Mass = mass,
Chamber = chamber,
System = system1,
Volume = chamber_vol,
Date_tested = Date_tested,
Date_analysed =Date_analysed,
Swim = Swim,
Salinity = salinity,
Temperature = as.numeric(unique(firesting2$temperature)),
Resting_kg = mean(apoly_rmr$rate.output*-1),
Resting = mean(apoly_rmr$rate.output*-1)*mass,
rsqrest =mean(apoly_rmr$rsq))
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 59 | 1 | CARL354 | CARL355 | Arlington reef | 281 | 0.0006711 | ch4 | Dell | 0.0486 | 2023-04-04 | 2024-06-26 | good/good | 36 | 27 | 177.4127 | 0.1190617 | 0.9658 |
Maximum oxygen consumption
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch4
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 0.94 1.89
## -----------------------------------------
Subset data
cycle1.start <- Cycle_1.mmr[1,1]
cycle1.end <- tail(Cycle_1.mmr, n=1)[1,1]
cycle1.start.row <- which(firesting2_mmr$TIME == cycle1.start); cycle1.start## Warning in which(firesting2_mmr$TIME == cycle1.start): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="## Warning in which(firesting2_mmr$TIME == cycle1.end): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="cycle1_data <- firesting2_mmr |>
subset_data(from = cycle1.start.row,
to = cycle1.end.row,
by = "row") ## subset_data: Multi-column dataset detected in input!
## subset_data is generally intended to subset data already passed through inspect(), or 2-column data frames where time and oxygen are in columns 1 and 2 respectively.
## Subsetting will proceed anyway based on this assumption, but please ensure you understand what you are doing.## inspect: Applying column default of 'time = 1'## inspect: Applying column default of 'oxygen = 2'## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch4
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 2 3 4 5 6 7 8 9 10 11 12 14 15 16 17 18 19 20 21 22
## Minimum and Maximum intervals in uneven Time data:
## [1] 0.94 1.41
## -----------------------------------------
Calculating MMR
## Warning: auto_rate: Multi-column dataset detected in input. Selecting first two columns by default.
## If these are not the intended data, inspect() or subset the data frame columns appropriately before running auto_rate()
##
## # summary.auto_rate # -------------------
##
## === Summary of Results by Highest Rate ===
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 392.6590 -0.03827657 0.9857498 NA 136 196 7763.41
## 2: NA 2 392.3415 -0.03823576 0.9856901 NA 137 197 7764.38
## 3: NA 3 392.0464 -0.03819806 0.9856065 NA 135 195 7762.40
## 4: NA 4 391.7585 -0.03816088 0.9855650 NA 138 198 7765.37
## 5: NA 5 391.3135 -0.03810404 0.9856579 NA 139 199 7766.37
## ---
## 237: NA 237 219.8386 -0.01593271 0.9462477 NA 5 65 7632.46
## 238: NA 238 219.4768 -0.01588562 0.9455518 NA 4 64 7631.48
## 239: NA 239 217.6095 -0.01564200 0.9423784 NA 3 63 7630.46
## 240: NA 240 215.8673 -0.01541460 0.9402137 NA 2 62 7629.47
## 241: NA 241 212.9021 -0.01502761 0.9336263 NA 1 61 7628.47
## endtime oxy endoxy rate
## 1: 7823.41 95.455 93.207 -0.03827657
## 2: 7824.38 95.462 93.198 -0.03823576
## 3: 7822.40 95.452 93.179 -0.03819806
## 4: 7825.37 95.429 93.177 -0.03816088
## 5: 7826.37 95.458 93.060 -0.03810404
## ---
## 237: 7692.46 98.189 97.293 -0.01593271
## 238: 7691.48 98.145 97.248 -0.01588562
## 239: 7690.46 98.175 97.259 -0.01564200
## 240: 7689.47 98.105 97.262 -0.01541460
## 241: 7688.47 98.176 97.315 -0.01502761
##
## Regressions : 241 | Results : 241 | Method : highest | Roll width : 60 | Roll type : time
## -----------------------------------------Adjusting
## adjust_rate: Rate adjustments applied using "mean" method.##
## # print.adjust_rate # -------------------
## NOTE: Consider the sign of the adjustment value when adjusting the rate.
##
## Adjustment was applied using the 'mean' method.
##
## Rank 1 of 241 adjusted rate(s):
## Rate : -0.03827657
## Adjustment : 0.0004003393
## Adjusted Rate : -0.03867691
##
## To see other results use 'pos' input.
## To see full results use summary().
## -----------------------------------------Converting units
mmr_adj2 <- mmr_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253)## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.selecting rates
mmr_final <- mmr_adj2 |>
select_rate(method = "rsq", n=c(0.93,1)) |>
select_rate(method = "highest", n=1) |>
plot(type="full") |>
summary(export=TRUE)## select_rate: Selecting rates with rsq values between 0.93 and 1...## ----- Selection complete. 0 rate(s) removed, 241 rate(s) remaining -----## select_rate: Selecting highest 1 *absolute* rate values...## ----- Selection complete. 240 rate(s) removed, 1 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 392.659 -0.03827657 0.9857498 NA 136 196 7763.41
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 7823.41 95.455 93.207 -0.03827657 0.0004003393 -0.03867691 -0.03867691
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.0486 0.0006711 NA 36 27 1.013253 -0.4392912
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -654.5838 NA mgO2/hr/kg -654.5838
## -----------------------------------------Results
results <- results |>
mutate(Max_kg = mmr_final$rate.output*-1,
Max = (mmr_final$rate.output*-1)*mass,
rsqmax =mmr_final$rsq,
AAS_kg = Max_kg - Resting_kg,
AAS = Max - Resting,
Notes=Notes,
True_resting="")
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest | Max_kg | Max | rsqmax | AAS_kg | AAS | Notes | True_resting |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 59 | 1 | CARL354 | CARL355 | Arlington reef | 281 | 0.0006711 | ch4 | Dell | 0.0486 | 2023-04-04 | 2024-06-26 | good/good | 36 | 27 | 177.4127 | 0.1190617 | 0.9658 | 654.5838 | 0.4392912 | 0.9857498 | 477.1711 | 0.3202295 | check max |
Exporting data
## Rows: 265 Columns: 25
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (9): Male, Female, Population, Chamber, System, Date_tested, Swim, Note...
## dbl (16): Clutch, Replicate, Tank, Mass, Volume, Date_analysed, Salinity, Te...
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.2
Enter specimen data
Replicate = 2
mass = 0.0007828
chamber = "ch3"
Swim = "good/good"
chamber_vol = chamber3_dell
system1 = "Dell"
Notes=""
##--- time of trail ---##
experiment_mmr_date <- "04 April 2023 02 37PM/Oxygen"
experiment_mmr_date2 <- "04 April 2023 02 37PM/All"
firesting_mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_mmr_date,"data raw/Firesting.txt"),
delim = "\t", escape_double = FALSE,
col_types = cols(`Time (HH:MM:SS)` = col_time(format = "%H:%M:%S"),
`Time (s)` = col_number(), Ch1...5 = col_number(),
Ch2...6 = col_number(), Ch3...7 = col_number(),
Ch4...8 = col_number()), trim_ws = TRUE,
skip = 19) ## New names:
## • `Ch1` -> `Ch1...5`
## • `Ch2` -> `Ch2...6`
## • `Ch3` -> `Ch3...7`
## • `Ch4` -> `Ch4...8`
## • `Ch 1` -> `Ch 1...9`
## • `Ch 2` -> `Ch 2...10`
## • `Ch 3` -> `Ch 3...11`
## • `Ch 4` -> `Ch 4...12`
## • `('C)` -> `('C)...15`
## • `('C)` -> `('C)...16`
## • `Ch 1` -> `Ch 1...18`
## • `Ch 2` -> `Ch 2...19`
## • `Ch 3` -> `Ch 3...20`
## • `Ch 4` -> `Ch 4...21`
## • `Ch1` -> `Ch1...22`
## • `Ch2` -> `Ch2...23`
## • `Ch3` -> `Ch3...24`
## • `Ch4` -> `Ch4...25`
## • `Ch1` -> `Ch1...26`
## • `Ch2` -> `Ch2...27`
## • `Ch3` -> `Ch3...28`
## • `Ch4` -> `Ch4...29`
## • `` -> `...31`## Warning: One or more parsing issues, call `problems()` on your data frame for details,
## e.g.:
## dat <- vroom(...)
## problems(dat)Cycle_1.mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_mmr_date2,"slopes/Cycle_1.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE) ## New names:
## • `` -> `...8`Background rates
Pre-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",preexperiment_date,"slopes"))
pre_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_pre1 <- pre_cycle1 %>% calc_rate.bg()##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.0001705656post-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",postexperiment_date,"slopes"))
post_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_post1 <- post_cycle1 %>% calc_rate.bg() ##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] 0.001973166Resting metabolic rate
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch3
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 0.94 3.86
## -----------------------------------------
#### subset data
Tstart.row=which(firesting2$TIME ==Cycle_1$Time[1], firesting$TIME)[1]# custom
Tstart.dTIME=as.numeric(firesting2[Tstart.row, "dTIME"])
Tend.row=which(firesting2$TIME ==tail(Cycle_last$Time, n=1), firesting$TIME)
Tend.dTIME=as.numeric(firesting2[Tend.row, "dTIME"])
apoly_insp <- firesting2 |>
subset_data(from=Tstart.dTIME,
to=Tend.dTIME,
by="time")
apoly_insp <- inspect(apoly_insp, time=1, oxygen=2)##
## # print.inspect # -----------------------
## dTIME ch3
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 7 8 9 10 11 12 14 16 18 19 20 21 22 23 24
## Minimum and Maximum intervals in uneven Time data:
## [1] 0.94 3.86
## -----------------------------------------
Extract rates
apoly_cr.int <- calc_rate.int(apoly_insp,
starts=(195+45+300),
wait=45,
measure=255,
by="time",
plot=TRUE) ##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from all replicates ...
## -----------------------------------------##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from selected replicates...
## To plot others modify 'pos' input.
## -----------------------------------------adjust rates for background
apoly_cr.int_adj <- adjust_rate(apoly_cr.int,
by = bg_pre,
by2 = bg_post,
time_by = Tstart.row,
time_by2 = Tend.row,
method = "linear")## Warning: adjust_rate: One or more of the timestamps for the rate(s) in 'x' do not lie between the timestamps for the 'by' and 'by2' background rates.
## Ensure this is correct. The adjustment value has been calculated regardless by extrapolating outside the background rates time window.## Warning: adjust_rate: background rates in 'by' and 'by2' differ in sign (i.e. one is +ve, one is -ve).
## Ensure this is correct. The 'linear' adjustment has been performed regardless.## adjust_rate: Rate adjustments applied using "linear" method.Converting units
apoly_cr.int_adj2 <- apoly_cr.int_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253) ## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.Plot curve
ggplot(as.data.frame(apoly_cr.int_adj2$summary), aes(x=row, y=rate.output*-1)) +
geom_point() +
stat_smooth(method = "lm", formula = y~poly(x, 2), color="red") +
theme_classic()
Rate filtering
apoly_rmr <- apoly_cr.int_adj2 |>
select_rate(method ="rsq", n=c(0.95,1)) |>
select_rate(method="lowest", n=6) |>
plot(type="full") |>
summary(export = TRUE)## select_rate: Selecting rates with rsq values between 0.95 and 1...## ----- Selection complete. 3 rate(s) removed, 17 rate(s) remaining -----## select_rate: Selecting lowest 6 *absolute* rate values...## ----- Selection complete. 11 rate(s) removed, 6 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: 7 1 258.2240 -0.01248817 0.971 NA 3220 3475 12796.69
## 2: 10 1 261.3960 -0.01129043 0.973 NA 4840 5095 14416.69
## 3: 12 1 273.0543 -0.01122669 0.984 NA 5919 6175 15496.19
## 4: 13 1 281.6668 -0.01141315 0.978 NA 6460 6715 16036.69
## 5: 16 1 281.6299 -0.01036707 0.977 NA 8080 8335 17656.71
## 6: 18 1 286.8860 -0.01003351 0.989 NA 9160 9415 18736.69
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 13051.71 98.096 94.847 -0.01248817 0.0005925848 -0.01308075 -0.01308075
## 2: 14671.69 98.331 95.416 -0.01129043 0.0009233929 -0.01221382 -0.01221382
## 3: 15751.81 98.848 96.050 -0.01122669 0.0011438942 -0.01237058 -0.01237058
## 4: 16291.69 98.478 95.390 -0.01141315 0.0012542031 -0.01266735 -0.01266735
## 5: 17911.79 98.273 95.718 -0.01036707 0.0015850255 -0.01195209 -0.01195209
## 6: 18991.70 98.632 96.210 -0.01003351 0.0018055543 -0.01183907 -0.01183907
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04977 0.0007828 NA 36 27 1.013253 -0.1521475
## 2: %Air sec 0.04977 0.0007828 NA 36 27 1.013253 -0.1420639
## 3: %Air sec 0.04977 0.0007828 NA 36 27 1.013253 -0.1438873
## 4: %Air sec 0.04977 0.0007828 NA 36 27 1.013253 -0.1473391
## 5: %Air sec 0.04977 0.0007828 NA 36 27 1.013253 -0.1390196
## 6: %Air sec 0.04977 0.0007828 NA 36 27 1.013253 -0.1377050
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -194.3632 NA mgO2/hr/kg -194.3632
## 2: -181.4817 NA mgO2/hr/kg -181.4817
## 3: -183.8110 NA mgO2/hr/kg -183.8110
## 4: -188.2206 NA mgO2/hr/kg -188.2206
## 5: -177.5927 NA mgO2/hr/kg -177.5927
## 6: -175.9133 NA mgO2/hr/kg -175.9133
## -----------------------------------------Results
results <- data.frame(Clutch = Clutch,
Replicate =Replicate,
Male=Male,
Female=Female,
Population = Population,
Tank = Tank,
Mass = mass,
Chamber = chamber,
System = system1,
Volume = chamber_vol,
Date_tested = Date_tested,
Date_analysed =Date_analysed,
Swim = Swim,
Salinity = salinity,
Temperature = as.numeric(unique(firesting2$temperature)),
Resting_kg = mean(apoly_rmr$rate.output*-1),
Resting = mean(apoly_rmr$rate.output*-1)*mass,
rsqrest =mean(apoly_rmr$rsq))
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 59 | 2 | CARL354 | CARL355 | Arlington reef | 281 | 0.0007828 | ch3 | Dell | 0.04977 | 2023-04-04 | 2024-06-26 | good/good | 36 | 27 | 185.0938 | 0.1448915 | 0.9766 |
Maximum oxygen consumption
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch3
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 0.94 1.89
## -----------------------------------------
Subset data
cycle1.start <- Cycle_1.mmr[1,1]
cycle1.end <- tail(Cycle_1.mmr, n=1)[1,1]
cycle1.start.row <- which(firesting2_mmr$TIME == cycle1.start); cycle1.start## Warning in which(firesting2_mmr$TIME == cycle1.start): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="## Warning in which(firesting2_mmr$TIME == cycle1.end): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="cycle1_data <- firesting2_mmr |>
subset_data(from = cycle1.start.row,
to = cycle1.end.row,
by = "row") ## subset_data: Multi-column dataset detected in input!
## subset_data is generally intended to subset data already passed through inspect(), or 2-column data frames where time and oxygen are in columns 1 and 2 respectively.
## Subsetting will proceed anyway based on this assumption, but please ensure you understand what you are doing.## inspect: Applying column default of 'time = 1'## inspect: Applying column default of 'oxygen = 2'## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch3
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
## Minimum and Maximum intervals in uneven Time data:
## [1] 0.94 1.18
## -----------------------------------------
Calculating MMR
## Warning: auto_rate: Multi-column dataset detected in input. Selecting first two columns by default.
## If these are not the intended data, inspect() or subset the data frame columns appropriately before running auto_rate()
##
## # summary.auto_rate # -------------------
##
## === Summary of Results by Highest Rate ===
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 546.6813 -0.05322153 0.9919990 NA 133 194 8438.56
## 2: NA 2 546.3530 -0.05318324 0.9918811 NA 132 193 8437.53
## 3: NA 3 546.1414 -0.05315729 0.9917891 NA 134 195 8439.55
## 4: NA 4 545.5037 -0.05308151 0.9915472 NA 135 196 8440.54
## 5: NA 5 545.1383 -0.05303785 0.9913980 NA 136 197 8441.54
## ---
## 236: NA 236 211.7532 -0.01351659 0.9662676 NA 47 108 8352.68
## 237: NA 237 211.2501 -0.01345595 0.9659388 NA 43 104 8348.71
## 238: NA 238 210.3320 -0.01334695 0.9667582 NA 46 107 8351.63
## 239: NA 239 210.3201 -0.01334508 0.9665977 NA 44 105 8349.67
## 240: NA 240 209.6702 -0.01326778 0.9679224 NA 45 106 8350.64
## endtime oxy endoxy rate
## 1: 8498.56 97.479 94.485 -0.05322153
## 2: 8497.53 97.482 94.501 -0.05318324
## 3: 8499.55 97.455 94.474 -0.05315729
## 4: 8500.54 97.409 94.428 -0.05308151
## 5: 8501.54 97.338 94.381 -0.05303785
## ---
## 236: 8412.68 98.834 97.970 -0.01351659
## 237: 8408.71 99.023 98.078 -0.01345595
## 238: 8411.63 98.906 97.984 -0.01334695
## 239: 8409.67 98.974 98.094 -0.01334508
## 240: 8410.64 98.978 98.040 -0.01326778
##
## Regressions : 240 | Results : 240 | Method : highest | Roll width : 60 | Roll type : time
## -----------------------------------------Adjusting
## adjust_rate: Rate adjustments applied using "mean" method.##
## # print.adjust_rate # -------------------
## NOTE: Consider the sign of the adjustment value when adjusting the rate.
##
## Adjustment was applied using the 'mean' method.
##
## Rank 1 of 240 adjusted rate(s):
## Rate : -0.05322153
## Adjustment : -0.0001705656
## Adjusted Rate : -0.05305096
##
## To see other results use 'pos' input.
## To see full results use summary().
## -----------------------------------------Converting units
mmr_adj2 <- mmr_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253)## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.selecting rates
mmr_final <- mmr_adj2 |>
select_rate(method = "rsq", n=c(0.93,1)) |>
select_rate(method = "highest", n=1) |>
plot(type="full") |>
summary(export=TRUE)## select_rate: Selecting rates with rsq values between 0.93 and 1...## ----- Selection complete. 11 rate(s) removed, 229 rate(s) remaining -----## select_rate: Selecting highest 1 *absolute* rate values...## ----- Selection complete. 228 rate(s) removed, 1 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 546.6813 -0.05322153 0.991999 NA 133 194 8438.56
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 8498.56 97.479 94.485 -0.05322153 -0.0001705656 -0.05305096 -0.05305096
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04977 0.0007828 NA 36 27 1.013253 -0.6170572
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -788.2693 NA mgO2/hr/kg -788.2693
## -----------------------------------------Results
results <- results |>
mutate(Max_kg = mmr_final$rate.output*-1,
Max = (mmr_final$rate.output*-1)*mass,
rsqmax =mmr_final$rsq,
AAS_kg = Max_kg - Resting_kg,
AAS = Max - Resting,
Notes=Notes,
True_resting="")
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest | Max_kg | Max | rsqmax | AAS_kg | AAS | Notes | True_resting |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 59 | 2 | CARL354 | CARL355 | Arlington reef | 281 | 0.0007828 | ch3 | Dell | 0.04977 | 2023-04-04 | 2024-06-26 | good/good | 36 | 27 | 185.0938 | 0.1448915 | 0.9766 | 788.2693 | 0.6170572 | 0.991999 | 603.1754 | 0.4721657 |
Exporting data
## Rows: 266 Columns: 25
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (9): Male, Female, Population, Chamber, System, Date_tested, Swim, Note...
## dbl (16): Clutch, Replicate, Tank, Mass, Volume, Date_analysed, Salinity, Te...
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.3
Enter specimen data
Replicate = 3
mass = 0.0006218
chamber = "ch2"
Swim = "good/good"
chamber_vol = chamber2_dell
system1 = "Dell"
Notes="check max"
##--- time of trail ---##
experiment_mmr_date <- "04 April 2023 02 48PM/Oxygen"
experiment_mmr_date2 <- "04 April 2023 02 48PM/All"
firesting_mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_mmr_date,"data raw/Firesting.txt"),
delim = "\t", escape_double = FALSE,
col_types = cols(`Time (HH:MM:SS)` = col_time(format = "%H:%M:%S"),
`Time (s)` = col_number(), Ch1...5 = col_number(),
Ch2...6 = col_number(), Ch3...7 = col_number(),
Ch4...8 = col_number()), trim_ws = TRUE,
skip = 19) ## New names:
## • `Ch1` -> `Ch1...5`
## • `Ch2` -> `Ch2...6`
## • `Ch3` -> `Ch3...7`
## • `Ch4` -> `Ch4...8`
## • `Ch 1` -> `Ch 1...9`
## • `Ch 2` -> `Ch 2...10`
## • `Ch 3` -> `Ch 3...11`
## • `Ch 4` -> `Ch 4...12`
## • `('C)` -> `('C)...15`
## • `('C)` -> `('C)...16`
## • `Ch 1` -> `Ch 1...18`
## • `Ch 2` -> `Ch 2...19`
## • `Ch 3` -> `Ch 3...20`
## • `Ch 4` -> `Ch 4...21`
## • `Ch1` -> `Ch1...22`
## • `Ch2` -> `Ch2...23`
## • `Ch3` -> `Ch3...24`
## • `Ch4` -> `Ch4...25`
## • `Ch1` -> `Ch1...26`
## • `Ch2` -> `Ch2...27`
## • `Ch3` -> `Ch3...28`
## • `Ch4` -> `Ch4...29`
## • `` -> `...31`## Warning: One or more parsing issues, call `problems()` on your data frame for details,
## e.g.:
## dat <- vroom(...)
## problems(dat)Cycle_1.mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_mmr_date2,"slopes/Cycle_1.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE) ## New names:
## • `` -> `...8`Background rates
Pre-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",preexperiment_date,"slopes"))
pre_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_pre1 <- pre_cycle1 %>% calc_rate.bg()##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] 0.0001990685post-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",postexperiment_date,"slopes"))
post_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_post1 <- post_cycle1 %>% calc_rate.bg() ##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] 0.0001371767Resting metabolic rate
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch2
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 0.94 3.86
## -----------------------------------------
#### subset data
Tstart.row=which(firesting2$TIME ==Cycle_1$Time[1], firesting$TIME)[1]# custom
Tstart.dTIME=as.numeric(firesting2[Tstart.row, "dTIME"])
Tend.row=which(firesting2$TIME ==tail(Cycle_last$Time, n=1), firesting$TIME)
Tend.dTIME=as.numeric(firesting2[Tend.row, "dTIME"])
apoly_insp <- firesting2 |>
subset_data(from=Tstart.dTIME,
to=Tend.dTIME,
by="time")
apoly_insp <- inspect(apoly_insp, time=1, oxygen=2)##
## # print.inspect # -----------------------
## dTIME ch2
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 7 8 9 10 11 12 14 16 18 19 20 21 22 23 24
## Minimum and Maximum intervals in uneven Time data:
## [1] 0.94 3.86
## -----------------------------------------
Extract rates
apoly_cr.int <- calc_rate.int(apoly_insp,
starts=(195+45+300),
wait=45,
measure=255,
by="time",
plot=TRUE) ##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from all replicates ...
## -----------------------------------------##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from selected replicates...
## To plot others modify 'pos' input.
## -----------------------------------------adjust rates for background
apoly_cr.int_adj <- adjust_rate(apoly_cr.int,
by = bg_pre,
by2 = bg_post,
time_by = Tstart.row,
time_by2 = Tend.row,
method = "linear")## Warning: adjust_rate: One or more of the timestamps for the rate(s) in 'x' do not lie between the timestamps for the 'by' and 'by2' background rates.
## Ensure this is correct. The adjustment value has been calculated regardless by extrapolating outside the background rates time window.## adjust_rate: Rate adjustments applied using "linear" method.Converting units
apoly_cr.int_adj2 <- apoly_cr.int_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253) ## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.Plot curve
ggplot(as.data.frame(apoly_cr.int_adj2$summary), aes(x=row, y=rate.output*-1)) +
geom_point() +
stat_smooth(method = "lm", formula = y~poly(x, 2), color="red") +
theme_classic()
Rate filtering
apoly_rmr <- apoly_cr.int_adj2 |>
select_rate(method ="rsq", n=c(0.95,1)) |>
select_rate(method="lowest", n=6) |>
plot(type="full") |>
summary(export = TRUE)## select_rate: Selecting rates with rsq values between 0.95 and 1...## ----- Selection complete. 0 rate(s) removed, 20 rate(s) remaining -----## select_rate: Selecting lowest 6 *absolute* rate values...## ----- Selection complete. 14 rate(s) removed, 6 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: 2 1 271.6974 -0.01707716 0.977 NA 586 814 10096.80
## 2: 3 1 293.8113 -0.01824522 0.967 NA 1059 1315 10636.47
## 3: 4 1 286.7919 -0.01675359 0.978 NA 1600 1855 11176.71
## 4: 5 1 315.8395 -0.01844303 0.982 NA 2140 2395 11716.70
## 5: 15 1 397.2490 -0.01740853 0.992 NA 7540 7795 17116.78
## 6: 20 1 478.1063 -0.01912128 0.996 NA 10240 10495 19816.69
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 10351.57 98.868 94.380 -0.01707716 0.0001929537 -0.01727012 -0.01727012
## 2: 10891.69 99.199 94.478 -0.01824522 0.0001897707 -0.01843499 -0.01843499
## 3: 11431.69 99.215 94.856 -0.01675359 0.0001865864 -0.01694018 -0.01694018
## 4: 11971.69 99.233 94.752 -0.01844303 0.0001834028 -0.01862643 -0.01862643
## 5: 17371.72 99.016 94.569 -0.01740853 0.0001515664 -0.01756009 -0.01756009
## 6: 20071.70 99.059 94.120 -0.01912128 0.0001356486 -0.01925693 -0.01925693
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04593 0.0006218 NA 36 27 1.013253 -0.1853772
## 2: %Air sec 0.04593 0.0006218 NA 36 27 1.013253 -0.1978809
## 3: %Air sec 0.04593 0.0006218 NA 36 27 1.013253 -0.1818356
## 4: %Air sec 0.04593 0.0006218 NA 36 27 1.013253 -0.1999358
## 5: %Air sec 0.04593 0.0006218 NA 36 27 1.013253 -0.1884897
## 6: %Air sec 0.04593 0.0006218 NA 36 27 1.013253 -0.2067035
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -298.1299 NA mgO2/hr/kg -298.1299
## 2: -318.2388 NA mgO2/hr/kg -318.2388
## 3: -292.4342 NA mgO2/hr/kg -292.4342
## 4: -321.5437 NA mgO2/hr/kg -321.5437
## 5: -303.1357 NA mgO2/hr/kg -303.1357
## 6: -332.4277 NA mgO2/hr/kg -332.4277
## -----------------------------------------Results
results <- data.frame(Clutch = Clutch,
Replicate =Replicate,
Male=Male,
Female=Female,
Population = Population,
Tank = Tank,
Mass = mass,
Chamber = chamber,
System = system1,
Volume = chamber_vol,
Date_tested = Date_tested,
Date_analysed =Date_analysed,
Swim = Swim,
Salinity = salinity,
Temperature = as.numeric(unique(firesting2$temperature)),
Resting_kg = mean(apoly_rmr$rate.output*-1),
Resting = mean(apoly_rmr$rate.output*-1)*mass,
rsqrest =mean(apoly_rmr$rsq))
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 59 | 3 | CARL354 | CARL355 | Arlington reef | 281 | 0.0006218 | ch2 | Dell | 0.04593 | 2023-04-04 | 2024-06-26 | good/good | 36 | 27 | 314.6951 | 0.1956774 | 0.9828 |
Maximum oxygen consumption
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch2
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 0.94 1.94
## -----------------------------------------
Subset data
cycle1.start <- Cycle_1.mmr[1,1]
cycle1.end <- tail(Cycle_1.mmr, n=1)[1,1]
cycle1.start.row <- which(firesting2_mmr$TIME == cycle1.start); cycle1.start## Warning in which(firesting2_mmr$TIME == cycle1.start): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="## Warning in which(firesting2_mmr$TIME == cycle1.end): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="cycle1_data <- firesting2_mmr |>
subset_data(from = cycle1.start.row,
to = cycle1.end.row,
by = "row") ## subset_data: Multi-column dataset detected in input!
## subset_data is generally intended to subset data already passed through inspect(), or 2-column data frames where time and oxygen are in columns 1 and 2 respectively.
## Subsetting will proceed anyway based on this assumption, but please ensure you understand what you are doing.## inspect: Applying column default of 'time = 1'## inspect: Applying column default of 'oxygen = 2'## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch2
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 10 11 12 13 14 15 17 18 19 20 21 22
## Minimum and Maximum intervals in uneven Time data:
## [1] 0.94 1.17
## -----------------------------------------
Calculating MMR
## Warning: auto_rate: Multi-column dataset detected in input. Selecting first two columns by default.
## If these are not the intended data, inspect() or subset the data frame columns appropriately before running auto_rate()
##
## # summary.auto_rate # -------------------
##
## === Summary of Results by Highest Rate ===
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 530.3106 -0.04765723 0.9930365 NA 240 301 9159.72
## 2: NA 2 528.6796 -0.04748015 0.9924111 NA 239 300 9158.77
## 3: NA 3 527.3545 -0.04733644 0.9917785 NA 238 299 9157.74
## 4: NA 4 526.4070 -0.04723378 0.9913638 NA 237 298 9156.78
## 5: NA 5 525.1746 -0.04710019 0.9907093 NA 236 297 9155.77
## ---
## 236: NA 236 200.6981 -0.01134219 0.9338388 NA 29 90 8948.58
## 237: NA 237 200.6218 -0.01133188 0.9352826 NA 25 86 8944.46
## 238: NA 238 199.8782 -0.01125047 0.9413354 NA 28 89 8947.44
## 239: NA 239 199.6623 -0.01122548 0.9410584 NA 26 87 8945.48
## 240: NA 240 199.0969 -0.01116303 0.9448566 NA 27 88 8946.47
## endtime oxy endoxy rate
## 1: 9219.72 93.637 90.984 -0.04765723
## 2: 9218.77 93.702 91.067 -0.04748015
## 3: 9217.74 93.704 91.114 -0.04733644
## 4: 9216.78 93.746 91.176 -0.04723378
## 5: 9215.77 93.733 91.164 -0.04710019
## ---
## 236: 9008.58 99.266 98.388 -0.01134219
## 237: 9004.46 99.454 98.546 -0.01133188
## 238: 9007.44 99.307 98.408 -0.01125047
## 239: 9005.48 99.433 98.465 -0.01122548
## 240: 9006.47 99.409 98.430 -0.01116303
##
## Regressions : 240 | Results : 240 | Method : highest | Roll width : 60 | Roll type : time
## -----------------------------------------Adjusting
## adjust_rate: Rate adjustments applied using "mean" method.##
## # print.adjust_rate # -------------------
## NOTE: Consider the sign of the adjustment value when adjusting the rate.
##
## Adjustment was applied using the 'mean' method.
##
## Rank 1 of 240 adjusted rate(s):
## Rate : -0.04765723
## Adjustment : 0.0001990685
## Adjusted Rate : -0.0478563
##
## To see other results use 'pos' input.
## To see full results use summary().
## -----------------------------------------Converting units
mmr_adj2 <- mmr_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253)## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.selecting rates
mmr_final <- mmr_adj2 |>
select_rate(method = "rsq", n=c(0.93,1)) |>
select_rate(method = "highest", n=1) |>
plot(type="full") |>
summary(export=TRUE)## select_rate: Selecting rates with rsq values between 0.93 and 1...## ----- Selection complete. 18 rate(s) removed, 222 rate(s) remaining -----## select_rate: Selecting highest 1 *absolute* rate values...## ----- Selection complete. 221 rate(s) removed, 1 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 530.3106 -0.04765723 0.9930365 NA 240 301 9159.72
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 9219.72 93.637 90.984 -0.04765723 0.0001990685 -0.0478563 -0.0478563
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04593 0.0006218 NA 36 27 1.013253 -0.5136887
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -826.1317 NA mgO2/hr/kg -826.1317
## -----------------------------------------Results
results <- results |>
mutate(Max_kg = mmr_final$rate.output*-1,
Max = (mmr_final$rate.output*-1)*mass,
rsqmax =mmr_final$rsq,
AAS_kg = Max_kg - Resting_kg,
AAS = Max - Resting,
Notes=Notes,
True_resting="")
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest | Max_kg | Max | rsqmax | AAS_kg | AAS | Notes | True_resting |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 59 | 3 | CARL354 | CARL355 | Arlington reef | 281 | 0.0006218 | ch2 | Dell | 0.04593 | 2023-04-04 | 2024-06-26 | good/good | 36 | 27 | 314.6951 | 0.1956774 | 0.9828 | 826.1317 | 0.5136887 | 0.9930365 | 511.4366 | 0.3180113 | check max |
Exporting data
## Rows: 267 Columns: 25
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (9): Male, Female, Population, Chamber, System, Date_tested, Swim, Note...
## dbl (16): Clutch, Replicate, Tank, Mass, Volume, Date_analysed, Salinity, Te...
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.4
Enter specimen data
Replicate = 4
mass = 0.0004583
chamber = "ch1"
Swim = "good/good"
chamber_vol = chamber1_dell
system1 = "Dell"
Notes=""
##--- time of trail ---##
experiment_mmr_date <- "04 April 2023 02 58PM/Oxygen"
experiment_mmr_date2 <- "04 April 2023 02 58PM/All"
firesting_mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_mmr_date,"data raw/Firesting.txt"),
delim = "\t", escape_double = FALSE,
col_types = cols(`Time (HH:MM:SS)` = col_time(format = "%H:%M:%S"),
`Time (s)` = col_number(), Ch1...5 = col_number(),
Ch2...6 = col_number(), Ch3...7 = col_number(),
Ch4...8 = col_number()), trim_ws = TRUE,
skip = 19) ## New names:
## • `Ch1` -> `Ch1...5`
## • `Ch2` -> `Ch2...6`
## • `Ch3` -> `Ch3...7`
## • `Ch4` -> `Ch4...8`
## • `Ch 1` -> `Ch 1...9`
## • `Ch 2` -> `Ch 2...10`
## • `Ch 3` -> `Ch 3...11`
## • `Ch 4` -> `Ch 4...12`
## • `('C)` -> `('C)...15`
## • `('C)` -> `('C)...16`
## • `Ch 1` -> `Ch 1...18`
## • `Ch 2` -> `Ch 2...19`
## • `Ch 3` -> `Ch 3...20`
## • `Ch 4` -> `Ch 4...21`
## • `Ch1` -> `Ch1...22`
## • `Ch2` -> `Ch2...23`
## • `Ch3` -> `Ch3...24`
## • `Ch4` -> `Ch4...25`
## • `Ch1` -> `Ch1...26`
## • `Ch2` -> `Ch2...27`
## • `Ch3` -> `Ch3...28`
## • `Ch4` -> `Ch4...29`
## • `` -> `...31`## Warning: One or more parsing issues, call `problems()` on your data frame for details,
## e.g.:
## dat <- vroom(...)
## problems(dat)Cycle_1.mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_mmr_date2,"slopes/Cycle_1.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE) ## New names:
## • `` -> `...8`Background rates
Pre-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",preexperiment_date,"slopes"))
pre_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_pre1 <- pre_cycle1 %>% calc_rate.bg()##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.001411008post-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",postexperiment_date,"slopes"))
post_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_post1 <- post_cycle1 %>% calc_rate.bg() ##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.001481708Resting metabolic rate
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch1
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 0.94 3.86
## -----------------------------------------
#### subset data
Tstart.row=which(firesting2$TIME ==Cycle_1$Time[1], firesting$TIME)[1]# custom
Tstart.dTIME=as.numeric(firesting2[Tstart.row, "dTIME"])
Tend.row=which(firesting2$TIME ==tail(Cycle_last$Time, n=1), firesting$TIME)
Tend.dTIME=as.numeric(firesting2[Tend.row, "dTIME"])
apoly_insp <- firesting2 |>
subset_data(from=Tstart.dTIME,
to=Tend.dTIME,
by="time")
apoly_insp <- inspect(apoly_insp, time=1, oxygen=2)##
## # print.inspect # -----------------------
## dTIME ch1
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 7 8 9 10 11 12 14 16 18 19 20 21 22 23 24
## Minimum and Maximum intervals in uneven Time data:
## [1] 0.94 3.86
## -----------------------------------------
Extract rates
apoly_cr.int <- calc_rate.int(apoly_insp,
starts=(195+45+300),
wait=45,
measure=255,
by="time",
plot=TRUE) ##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from all replicates ...
## -----------------------------------------##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from selected replicates...
## To plot others modify 'pos' input.
## -----------------------------------------adjust rates for background
apoly_cr.int_adj <- adjust_rate(apoly_cr.int,
by = bg_pre,
by2 = bg_post,
time_by = Tstart.row,
time_by2 = Tend.row,
method = "linear")## Warning: adjust_rate: One or more of the timestamps for the rate(s) in 'x' do not lie between the timestamps for the 'by' and 'by2' background rates.
## Ensure this is correct. The adjustment value has been calculated regardless by extrapolating outside the background rates time window.## adjust_rate: Rate adjustments applied using "linear" method.Converting units
apoly_cr.int_adj2 <- apoly_cr.int_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253) ## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.Plot curve
ggplot(as.data.frame(apoly_cr.int_adj2$summary), aes(x=row, y=rate.output*-1)) +
geom_point() +
stat_smooth(method = "lm", formula = y~poly(x, 2), color="red") +
theme_classic()
Rate filtering
apoly_rmr <- apoly_cr.int_adj2 |>
select_rate(method ="rsq", n=c(0.95,1)) |>
select_rate(method="lowest", n=6) |>
plot(type="full") |>
summary(export = TRUE)## select_rate: Selecting rates with rsq values between 0.95 and 1...## ----- Selection complete. 0 rate(s) removed, 20 rate(s) remaining -----## select_rate: Selecting lowest 6 *absolute* rate values...## ----- Selection complete. 14 rate(s) removed, 6 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: 4 1 299.1518 -0.01788401 0.992 NA 1600 1855 11176.71
## 2: 12 1 387.3032 -0.01854843 0.993 NA 5919 6175 15496.19
## 3: 13 1 394.0143 -0.01834723 0.996 NA 6460 6715 16036.69
## 4: 14 1 405.2342 -0.01842554 0.987 NA 7000 7255 16576.69
## 5: 16 1 427.1591 -0.01854639 0.992 NA 8080 8335 17656.71
## 6: 18 1 415.1555 -0.01684735 0.974 NA 9160 9415 18736.69
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 11431.69 99.344 94.642 -0.01788401 -0.001425267 -0.01645875 -0.01645875
## 2: 15751.81 99.582 95.007 -0.01854843 -0.001454359 -0.01709407 -0.01709407
## 3: 16291.69 99.838 95.132 -0.01834723 -0.001457997 -0.01688923 -0.01688923
## 4: 16831.69 99.793 94.916 -0.01842554 -0.001461634 -0.01696391 -0.01696391
## 5: 17911.79 99.455 95.042 -0.01854639 -0.001468907 -0.01707749 -0.01707749
## 6: 18991.70 98.836 94.984 -0.01684735 -0.001476180 -0.01537117 -0.01537117
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.0465 0.0004583 NA 36 27 1.013253 -0.1788604
## 2: %Air sec 0.0465 0.0004583 NA 36 27 1.013253 -0.1857646
## 3: %Air sec 0.0465 0.0004583 NA 36 27 1.013253 -0.1835385
## 4: %Air sec 0.0465 0.0004583 NA 36 27 1.013253 -0.1843501
## 5: %Air sec 0.0465 0.0004583 NA 36 27 1.013253 -0.1855844
## 6: %Air sec 0.0465 0.0004583 NA 36 27 1.013253 -0.1670415
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -390.2693 NA mgO2/hr/kg -390.2693
## 2: -405.3341 NA mgO2/hr/kg -405.3341
## 3: -400.4769 NA mgO2/hr/kg -400.4769
## 4: -402.2476 NA mgO2/hr/kg -402.2476
## 5: -404.9408 NA mgO2/hr/kg -404.9408
## 6: -364.4807 NA mgO2/hr/kg -364.4807
## -----------------------------------------Results
results <- data.frame(Clutch = Clutch,
Replicate =Replicate,
Male=Male,
Female=Female,
Population = Population,
Tank = Tank,
Mass = mass,
Chamber = chamber,
System = system1,
Volume = chamber_vol,
Date_tested = Date_tested,
Date_analysed =Date_analysed,
Swim = Swim,
Salinity = salinity,
Temperature = as.numeric(unique(firesting2$temperature)),
Resting_kg = mean(apoly_rmr$rate.output*-1),
Resting = mean(apoly_rmr$rate.output*-1)*mass,
rsqrest =mean(apoly_rmr$rsq))
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 59 | 4 | CARL354 | CARL355 | Arlington reef | 281 | 0.0004583 | ch1 | Dell | 0.0465 | 2023-04-04 | 2024-06-26 | good/good | 36 | 27 | 400.6537 | 0.1836196 | 0.992 |
Maximum oxygen consumption
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch1
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 0.94 3.86
## -----------------------------------------
Subset data
cycle1.start <- Cycle_1.mmr[1,1]
cycle1.end <- tail(Cycle_1.mmr, n=1)[1,1]
cycle1.start.row <- which(firesting2_mmr$TIME == cycle1.start); cycle1.start## Warning in which(firesting2_mmr$TIME == cycle1.start): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="## Warning in which(firesting2_mmr$TIME == cycle1.end): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="cycle1_data <- firesting2_mmr |>
subset_data(from = cycle1.start.row[1], # custom
to = cycle1.end.row,
by = "row") ## subset_data: Multi-column dataset detected in input!
## subset_data is generally intended to subset data already passed through inspect(), or 2-column data frames where time and oxygen are in columns 1 and 2 respectively.
## Subsetting will proceed anyway based on this assumption, but please ensure you understand what you are doing.## inspect: Applying column default of 'time = 1'## inspect: Applying column default of 'oxygen = 2'## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch1
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 7 8 9 10 11 12 14 16 18 19 20 21 22 23 24
## Minimum and Maximum intervals in uneven Time data:
## [1] 0.94 1.21
## -----------------------------------------
Calculating MMR
## Warning: auto_rate: Multi-column dataset detected in input. Selecting first two columns by default.
## If these are not the intended data, inspect() or subset the data frame columns appropriately before running auto_rate()
##
## # summary.auto_rate # -------------------
##
## === Summary of Results by Highest Rate ===
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 541.5865 -0.04604043 0.9938490 NA 166 227 9677.73
## 2: NA 2 541.2575 -0.04600664 0.9938086 NA 165 226 9676.72
## 3: NA 3 539.1105 -0.04578563 0.9933582 NA 164 225 9675.72
## 4: NA 4 538.7102 -0.04574397 0.9932651 NA 167 228 9678.72
## 5: NA 5 537.7837 -0.04564897 0.9932545 NA 163 224 9674.72
## ---
## 237: NA 237 218.5470 -0.01259520 0.7912842 NA 62 123 9573.72
## 238: NA 238 218.2002 -0.01256003 0.7922028 NA 64 125 9575.72
## 239: NA 239 218.1526 -0.01255394 0.7908137 NA 61 122 9572.72
## 240: NA 240 217.5254 -0.01249019 0.7925169 NA 65 126 9576.72
## 241: NA 241 217.5277 -0.01248876 0.7895739 NA 60 121 9571.72
## endtime oxy endoxy rate
## 1: 9737.73 96.024 93.277 -0.04604043
## 2: 9736.72 95.980 93.295 -0.04600664
## 3: 9735.72 96.087 93.329 -0.04578563
## 4: 9738.72 96.053 93.364 -0.04574397
## 5: 9734.72 96.097 93.387 -0.04564897
## ---
## 237: 9633.72 98.022 97.139 -0.01259520
## 238: 9635.72 98.128 97.063 -0.01256003
## 239: 9632.72 97.998 97.176 -0.01255394
## 240: 9636.72 98.065 97.066 -0.01249019
## 241: 9631.72 98.067 97.260 -0.01248876
##
## Regressions : 241 | Results : 241 | Method : highest | Roll width : 60 | Roll type : time
## -----------------------------------------Adjusting
## adjust_rate: Rate adjustments applied using "mean" method.##
## # print.adjust_rate # -------------------
## NOTE: Consider the sign of the adjustment value when adjusting the rate.
##
## Adjustment was applied using the 'mean' method.
##
## Rank 1 of 241 adjusted rate(s):
## Rate : -0.04604043
## Adjustment : -0.001411008
## Adjusted Rate : -0.04462942
##
## To see other results use 'pos' input.
## To see full results use summary().
## -----------------------------------------Converting units
mmr_adj2 <- mmr_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253)## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.selecting rates
mmr_final <- mmr_adj2 |>
select_rate(method = "rsq", n=c(0.93,1)) |>
select_rate(method = "highest", n=1) |>
plot(type="full") |>
summary(export=TRUE)## select_rate: Selecting rates with rsq values between 0.93 and 1...## ----- Selection complete. 56 rate(s) removed, 185 rate(s) remaining -----## select_rate: Selecting highest 1 *absolute* rate values...## ----- Selection complete. 184 rate(s) removed, 1 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 541.5865 -0.04604043 0.993849 NA 166 227 9677.73
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 9737.73 96.024 93.277 -0.04604043 -0.001411008 -0.04462942 -0.04462942
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.0465 0.0004583 NA 36 27 1.013253 -0.4849966
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -1058.251 NA mgO2/hr/kg -1058.251
## -----------------------------------------Results
results <- results |>
mutate(Max_kg = mmr_final$rate.output*-1,
Max = (mmr_final$rate.output*-1)*mass,
rsqmax =mmr_final$rsq,
AAS_kg = Max_kg - Resting_kg,
AAS = Max - Resting,
Notes=Notes,
True_resting="")
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest | Max_kg | Max | rsqmax | AAS_kg | AAS | Notes | True_resting |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 59 | 4 | CARL354 | CARL355 | Arlington reef | 281 | 0.0004583 | ch1 | Dell | 0.0465 | 2023-04-04 | 2024-06-26 | good/good | 36 | 27 | 400.6537 | 0.1836196 | 0.992 | 1058.251 | 0.4849966 | 0.993849 | 657.5975 | 0.301377 |
Exporting data
## Rows: 268 Columns: 25
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (9): Male, Female, Population, Chamber, System, Date_tested, Swim, Note...
## dbl (16): Clutch, Replicate, Tank, Mass, Volume, Date_analysed, Salinity, Te...
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.5
Enter specimen data
Replicate = 5
mass = 0.0005481
chamber = "ch4"
Swim = "good/good"
chamber_vol = chamber4_asus
system1 = "Asus"
Notes="check max; may not be reliable"
##--- time of trail ---##
experiment_mmr_date_asus <- "04 April 2023 02 00PM/Oxygen"
experiment_mmr_date2_asus <- "04 April 2023 02 00PM/All"
firesting_mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_mmr_date_asus,"data raw/Firesting.txt"),
delim = "\t", escape_double = FALSE,
col_types = cols(`Time (HH:MM:SS)` = col_time(format = "%H:%M:%S"),
`Time (s)` = col_number(), Ch1...5 = col_number(),
Ch2...6 = col_number(), Ch3...7 = col_number(),
Ch4...8 = col_number()), trim_ws = TRUE,
skip = 19) ## New names:
## • `Ch1` -> `Ch1...5`
## • `Ch2` -> `Ch2...6`
## • `Ch3` -> `Ch3...7`
## • `Ch4` -> `Ch4...8`
## • `Ch 1` -> `Ch 1...9`
## • `Ch 2` -> `Ch 2...10`
## • `Ch 3` -> `Ch 3...11`
## • `Ch 4` -> `Ch 4...12`
## • `('C)` -> `('C)...15`
## • `('C)` -> `('C)...16`
## • `Ch 1` -> `Ch 1...18`
## • `Ch 2` -> `Ch 2...19`
## • `Ch 3` -> `Ch 3...20`
## • `Ch 4` -> `Ch 4...21`
## • `Ch1` -> `Ch1...22`
## • `Ch2` -> `Ch2...23`
## • `Ch3` -> `Ch3...24`
## • `Ch4` -> `Ch4...25`
## • `Ch1` -> `Ch1...26`
## • `Ch2` -> `Ch2...27`
## • `Ch3` -> `Ch3...28`
## • `Ch4` -> `Ch4...29`
## • `` -> `...31`## Warning: One or more parsing issues, call `problems()` on your data frame for details,
## e.g.:
## dat <- vroom(...)
## problems(dat)Cycle_1.mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_mmr_date2_asus,"slopes/Cycle_1.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE) ## New names:
## • `` -> `...8`Background rates
Pre-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",preexperiment_date_asus,"slopes"))
pre_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_pre1 <- pre_cycle1 %>% calc_rate.bg()##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.0004067376post-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",postexperiment_date_asus,"slopes"))
post_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_post1 <- post_cycle1 %>% calc_rate.bg() ##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.0009798081Resting metabolic rate
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch4
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.06 4.06
## -----------------------------------------
#### subset data
Tstart.row=which(firesting2_asus$TIME ==Cycle_1_asus$Time[1], firesting_asus$TIME)
Tstart.dTIME=as.numeric(firesting2_asus[Tstart.row, "dTIME"])
Tend.row=which(firesting2_asus$TIME ==tail(Cycle_last_asus$Time, n=1), firesting_asus$TIME)
Tend.dTIME=as.numeric(firesting2_asus[Tend.row, "dTIME"])
apoly_insp <- firesting2_asus |>
subset_data(from=Tstart.dTIME,
to=Tend.dTIME,
by="time")
apoly_insp <- inspect(apoly_insp, time=1, oxygen=2)##
## # print.inspect # -----------------------
## dTIME ch4
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 3 6 7 8 9 10 11 13 14 16 17 19 20 22 24 25 26 28 29
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.06 3.52
## -----------------------------------------
Extract rates
apoly_cr.int <- calc_rate.int(apoly_insp,
starts=(195+45+300),
wait=45,
measure=255,
by="time",
plot=TRUE) ##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from all replicates ...
## -----------------------------------------adjust rates for background
apoly_cr.int_adj <- adjust_rate(apoly_cr.int,
by = bg_pre,
by2 = bg_post,
time_by = Tstart.row,
time_by2 = Tend.row,
method = "linear")## Warning: adjust_rate: One or more of the timestamps for the rate(s) in 'x' do not lie between the timestamps for the 'by' and 'by2' background rates.
## Ensure this is correct. The adjustment value has been calculated regardless by extrapolating outside the background rates time window.## adjust_rate: Rate adjustments applied using "linear" method.Converting units
apoly_cr.int_adj2 <- apoly_cr.int_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2_asus$temperature)),
P = 1.013253) ## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.Plot curve
ggplot(as.data.frame(apoly_cr.int_adj2$summary), aes(x=row, y=rate.output*-1)) +
geom_point() +
stat_smooth(method = "lm", formula = y~poly(x, 2), color="red") +
theme_classic()
Rate filtering
apoly_rmr <- apoly_cr.int_adj2 |>
select_rate(method ="rsq", n=c(0.95,1)) |>
select_rate(method="lowest", n=6) |>
plot(type="full") |>
summary(export = TRUE)## select_rate: Selecting rates with rsq values between 0.95 and 1...## ----- Selection complete. 5 rate(s) removed, 15 rate(s) remaining -----## select_rate: Selecting lowest 6 *absolute* rate values...## ----- Selection complete. 9 rate(s) removed, 6 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: 4 1 226.8297 -0.01377932 0.974 NA 1552 1789 9369.49
## 2: 8 1 232.1376 -0.01155904 0.963 NA 3563 3801 11529.23
## 3: 9 1 264.0105 -0.01364832 0.965 NA 4066 4304 12068.62
## 4: 12 1 294.5756 -0.01426552 0.974 NA 5576 5814 13689.18
## 5: 13 1 301.0682 -0.01417561 0.966 NA 6080 6318 14229.55
## 6: 16 1 309.5398 -0.01327192 0.982 NA 7590 7828 15848.89
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 9623.96 97.453 93.675 -0.01377932 -0.0005425484 -0.01323677 -0.01323677
## 2: 11784.39 98.607 95.585 -0.01155904 -0.0006683621 -0.01089068 -0.01089068
## 3: 12323.78 98.772 95.423 -0.01364832 -0.0006997788 -0.01294855 -0.01294855
## 4: 13944.21 98.929 95.348 -0.01426552 -0.0007941642 -0.01347135 -0.01347135
## 5: 14484.59 98.954 95.332 -0.01417561 -0.0008256382 -0.01334997 -0.01334997
## 6: 16103.92 98.978 95.575 -0.01327192 -0.0009199560 -0.01235196 -0.01235196
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.05068 0.0005481 NA 36 27 1.013253 -0.1567773
## 2: %Air sec 0.05068 0.0005481 NA 36 27 1.013253 -0.1289900
## 3: %Air sec 0.05068 0.0005481 NA 36 27 1.013253 -0.1533635
## 4: %Air sec 0.05068 0.0005481 NA 36 27 1.013253 -0.1595557
## 5: %Air sec 0.05068 0.0005481 NA 36 27 1.013253 -0.1581181
## 6: %Air sec 0.05068 0.0005481 NA 36 27 1.013253 -0.1462976
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -286.0378 NA mgO2/hr/kg -286.0378
## 2: -235.3403 NA mgO2/hr/kg -235.3403
## 3: -279.8094 NA mgO2/hr/kg -279.8094
## 4: -291.1069 NA mgO2/hr/kg -291.1069
## 5: -288.4840 NA mgO2/hr/kg -288.4840
## 6: -266.9177 NA mgO2/hr/kg -266.9177
## -----------------------------------------Results
results <- data.frame(Clutch = Clutch,
Replicate =Replicate,
Male=Male,
Female=Female,
Population = Population,
Tank = Tank,
Mass = mass,
Chamber = chamber,
System = system1,
Volume = chamber_vol,
Date_tested = Date_tested,
Date_analysed =Date_analysed,
Swim = Swim,
Salinity = salinity,
Temperature = as.numeric(unique(firesting2$temperature)),
Resting_kg = mean(apoly_rmr$rate.output*-1),
Resting = mean(apoly_rmr$rate.output*-1)*mass,
rsqrest =mean(apoly_rmr$rsq))
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 59 | 5 | CARL354 | CARL355 | Arlington reef | 281 | 0.0005481 | ch4 | Asus | 0.05068 | 2023-04-04 | 2024-06-26 | good/good | 36 | 27 | 282.4711 | 0.1548224 | 0.9722 |
Maximum oxygen consumption
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch4
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.06 4.06
## -----------------------------------------
Subset data
cycle1.start <- Cycle_1.mmr[1,1]
cycle1.end <- tail(Cycle_1.mmr, n=1)[1,1]
cycle1.start.row <- which(firesting2_mmr$TIME == cycle1.start); cycle1.start## Warning in which(firesting2_mmr$TIME == cycle1.start): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="## Warning in which(firesting2_mmr$TIME == cycle1.end): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="cycle1_data <- firesting2_mmr |>
subset_data(from = cycle1.start.row,
to = cycle1.end.row,
by = "row") ## subset_data: Multi-column dataset detected in input!
## subset_data is generally intended to subset data already passed through inspect(), or 2-column data frames where time and oxygen are in columns 1 and 2 respectively.
## Subsetting will proceed anyway based on this assumption, but please ensure you understand what you are doing.## inspect: Applying column default of 'time = 1'## inspect: Applying column default of 'oxygen = 2'## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch4
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 3 4 7 10 11 12 17 18 21 22 24 25 28 30 32 33 35 38 39 40
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.06 1.11
## -----------------------------------------
Calculating MMR
## Warning: auto_rate: Multi-column dataset detected in input. Selecting first two columns by default.
## If these are not the intended data, inspect() or subset the data frame columns appropriately before running auto_rate()
##
## # summary.auto_rate # -------------------
##
## === Summary of Results by Highest Rate ===
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 290.2122 -0.02697120 0.9868985 NA 32 88 7112.13
## 2: NA 2 289.9598 -0.02693626 0.9866100 NA 31 87 7111.05
## 3: NA 3 289.8529 -0.02691989 0.9865259 NA 38 94 7118.56
## 4: NA 4 289.7338 -0.02690348 0.9864363 NA 37 93 7117.49
## 5: NA 5 289.6032 -0.02688576 0.9865429 NA 33 89 7113.20
## ---
## 220: NA 220 180.0617 -0.01172112 0.9762301 NA 219 275 7312.91
## 221: NA 221 179.9648 -0.01170821 0.9738902 NA 223 279 7317.22
## 222: NA 222 179.7182 -0.01167450 0.9741858 NA 222 278 7316.14
## 223: NA 223 179.3646 -0.01162608 0.9748228 NA 220 276 7313.98
## 224: NA 224 179.3601 -0.01162558 0.9748231 NA 221 277 7315.05
## endtime oxy endoxy rate
## 1: 7172.13 98.299 96.842 -0.02697120
## 2: 7171.05 98.327 96.806 -0.02693626
## 3: 7178.56 98.169 96.651 -0.02691989
## 4: 7177.49 98.215 96.618 -0.02690348
## 5: 7173.20 98.351 96.804 -0.02688576
## ---
## 220: 7372.91 94.345 93.732 -0.01172112
## 221: 7377.22 94.314 93.554 -0.01170821
## 222: 7376.14 94.331 93.556 -0.01167450
## 223: 7373.98 94.338 93.682 -0.01162608
## 224: 7375.05 94.344 93.598 -0.01162558
##
## Regressions : 224 | Results : 224 | Method : highest | Roll width : 60 | Roll type : time
## -----------------------------------------Adjusting
## adjust_rate: Rate adjustments applied using "mean" method.##
## # print.adjust_rate # -------------------
## NOTE: Consider the sign of the adjustment value when adjusting the rate.
##
## Adjustment was applied using the 'mean' method.
##
## Rank 1 of 224 adjusted rate(s):
## Rate : -0.0269712
## Adjustment : -0.0004067376
## Adjusted Rate : -0.02656446
##
## To see other results use 'pos' input.
## To see full results use summary().
## -----------------------------------------Converting units
mmr_adj2 <- mmr_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253)## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.selecting rates
mmr_final <- mmr_adj2 |>
select_rate(method = "rsq", n=c(0.93,1)) |>
select_rate(method = "highest", n=1) |>
plot(type="full") |>
summary(export=TRUE)## select_rate: Selecting rates with rsq values between 0.93 and 1...## ----- Selection complete. 0 rate(s) removed, 224 rate(s) remaining -----## select_rate: Selecting highest 1 *absolute* rate values...## ----- Selection complete. 223 rate(s) removed, 1 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 290.2122 -0.0269712 0.9868985 NA 32 88 7112.13
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 7172.13 98.299 96.842 -0.0269712 -0.0004067376 -0.02656446 -0.02656446
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.05068 0.0005481 NA 36 27 1.013253 -0.3146314
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -574.0402 NA mgO2/hr/kg -574.0402
## -----------------------------------------Results
results <- results |>
mutate(Max_kg = mmr_final$rate.output*-1,
Max = (mmr_final$rate.output*-1)*mass,
rsqmax =mmr_final$rsq,
AAS_kg = Max_kg - Resting_kg,
AAS = Max - Resting,
Notes=Notes,
True_resting="")
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest | Max_kg | Max | rsqmax | AAS_kg | AAS | Notes | True_resting |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 59 | 5 | CARL354 | CARL355 | Arlington reef | 281 | 0.0005481 | ch4 | Asus | 0.05068 | 2023-04-04 | 2024-06-26 | good/good | 36 | 27 | 282.4711 | 0.1548224 | 0.9722 | 574.0402 | 0.3146314 | 0.9868985 | 291.569 | 0.159809 | check max; may not be reliable | |
| ### Expor | ting data |
## Rows: 269 Columns: 25
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (9): Male, Female, Population, Chamber, System, Date_tested, Swim, Note...
## dbl (16): Clutch, Replicate, Tank, Mass, Volume, Date_analysed, Salinity, Te...
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.6
Enter specimen data
Replicate = 6
mass = 0.0006868
chamber = "ch3"
Swim = "good/good"
chamber_vol = chamber3_asus
system1 = "Asus"
Notes=""
##--- time of trail ---##
experiment_mmr_date_asus <- "04 April 2023 01 11PM/Oxygen"
experiment_mmr_date2_asus <- "04 April 2023 01 11PM/All"
firesting_mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_mmr_date_asus,"data raw/Firesting.txt"),
delim = "\t", escape_double = FALSE,
col_types = cols(`Time (HH:MM:SS)` = col_time(format = "%H:%M:%S"),
`Time (s)` = col_number(), Ch1...5 = col_number(),
Ch2...6 = col_number(), Ch3...7 = col_number(),
Ch4...8 = col_number()), trim_ws = TRUE,
skip = 19) ## New names:
## • `Ch1` -> `Ch1...5`
## • `Ch2` -> `Ch2...6`
## • `Ch3` -> `Ch3...7`
## • `Ch4` -> `Ch4...8`
## • `Ch 1` -> `Ch 1...9`
## • `Ch 2` -> `Ch 2...10`
## • `Ch 3` -> `Ch 3...11`
## • `Ch 4` -> `Ch 4...12`
## • `('C)` -> `('C)...15`
## • `('C)` -> `('C)...16`
## • `Ch 1` -> `Ch 1...18`
## • `Ch 2` -> `Ch 2...19`
## • `Ch 3` -> `Ch 3...20`
## • `Ch 4` -> `Ch 4...21`
## • `Ch1` -> `Ch1...22`
## • `Ch2` -> `Ch2...23`
## • `Ch3` -> `Ch3...24`
## • `Ch4` -> `Ch4...25`
## • `Ch1` -> `Ch1...26`
## • `Ch2` -> `Ch2...27`
## • `Ch3` -> `Ch3...28`
## • `Ch4` -> `Ch4...29`
## • `` -> `...31`## Warning: One or more parsing issues, call `problems()` on your data frame for details,
## e.g.:
## dat <- vroom(...)
## problems(dat)Cycle_1.mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_mmr_date2_asus,"slopes/Cycle_1.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE) ## New names:
## • `` -> `...8`Background rates
Pre-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",preexperiment_date_asus,"slopes"))
pre_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_pre1 <- pre_cycle1 %>% calc_rate.bg()##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.0008082058post-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",postexperiment_date_asus,"slopes"))
post_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_post1 <- post_cycle1 %>% calc_rate.bg() ##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.001619877Resting metabolic rate
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch3
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.06 4.06
## -----------------------------------------
#### subset data
Tstart.row=which(firesting2_asus$TIME ==Cycle_1_asus$Time[1], firesting_asus$TIME)
Tstart.dTIME=as.numeric(firesting2_asus[Tstart.row, "dTIME"])
Tend.row=which(firesting2_asus$TIME ==tail(Cycle_last_asus$Time, n=1), firesting_asus$TIME)
Tend.dTIME=as.numeric(firesting2_asus[Tend.row, "dTIME"])
apoly_insp <- firesting2_asus |>
subset_data(from=Tstart.dTIME,
to=Tend.dTIME,
by="time")
apoly_insp <- inspect(apoly_insp, time=1, oxygen=2)##
## # print.inspect # -----------------------
## dTIME ch3
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 3 6 7 8 9 10 11 13 14 16 17 19 20 22 24 25 26 28 29
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.06 3.52
## -----------------------------------------
Extract rates
apoly_cr.int <- calc_rate.int(apoly_insp,
starts=(195+45+300),
wait=45,
measure=245,
by="time",
plot=TRUE) ##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from all replicates ...
## -----------------------------------------adjust rates for background
apoly_cr.int_adj <- adjust_rate(apoly_cr.int,
by = bg_pre,
by2 = bg_post,
time_by = Tstart.row,
time_by2 = Tend.row,
method = "linear")## Warning: adjust_rate: One or more of the timestamps for the rate(s) in 'x' do not lie between the timestamps for the 'by' and 'by2' background rates.
## Ensure this is correct. The adjustment value has been calculated regardless by extrapolating outside the background rates time window.## adjust_rate: Rate adjustments applied using "linear" method.Converting units
apoly_cr.int_adj2 <- apoly_cr.int_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253) ## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.Plot curve
ggplot(as.data.frame(apoly_cr.int_adj2$summary), aes(x=row, y=rate.output*-1)) +
geom_point() +
stat_smooth(method = "lm", formula = y~poly(x, 2), color="red") +
theme_classic()
Rate filtering
apoly_rmr <- apoly_cr.int_adj2 |>
select_rate(method ="rsq", n=c(0.95,1)) |>
select_rate(method="lowest", n=6) |>
plot(type="full") |>
summary(export = TRUE)## select_rate: Selecting rates with rsq values between 0.95 and 1...## ----- Selection complete. 0 rate(s) removed, 20 rate(s) remaining -----## select_rate: Selecting lowest 6 *absolute* rate values...## ----- Selection complete. 14 rate(s) removed, 6 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: 11 1 420.7884 -0.02441069 0.998 NA 5072 5301 13148.89
## 2: 13 1 450.6608 -0.02465678 0.999 NA 6080 6308 14229.55
## 3: 16 1 496.3438 -0.02502560 0.999 NA 7590 7819 15848.89
## 4: 17 1 514.4285 -0.02531611 0.999 NA 8094 8322 16389.30
## 5: 19 1 539.7029 -0.02519277 0.998 NA 9099 9327 17469.40
## 6: 20 1 527.5131 -0.02375634 0.999 NA 9602 9831 18008.80
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 13394.28 99.668 93.626 -0.02441069 -0.001311970 -0.02309872 -0.02309872
## 2: 14473.88 99.795 93.734 -0.02465678 -0.001401076 -0.02325570 -0.02325570
## 3: 16094.30 99.686 93.498 -0.02502560 -0.001534708 -0.02349089 -0.02349089
## 4: 16633.72 99.535 93.273 -0.02531611 -0.001579249 -0.02373686 -0.02373686
## 5: 17713.85 99.513 93.313 -0.02519277 -0.001668353 -0.02352442 -0.02352442
## 6: 18254.29 99.713 93.854 -0.02375634 -0.001712894 -0.02204344 -0.02204344
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04873 0.0006868 NA 36 27 1.013253 -0.2630564
## 2: %Air sec 0.04873 0.0006868 NA 36 27 1.013253 -0.2648441
## 3: %Air sec 0.04873 0.0006868 NA 36 27 1.013253 -0.2675225
## 4: %Air sec 0.04873 0.0006868 NA 36 27 1.013253 -0.2703237
## 5: %Air sec 0.04873 0.0006868 NA 36 27 1.013253 -0.2679044
## 6: %Air sec 0.04873 0.0006868 NA 36 27 1.013253 -0.2510385
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -383.0174 NA mgO2/hr/kg -383.0174
## 2: -385.6204 NA mgO2/hr/kg -385.6204
## 3: -389.5202 NA mgO2/hr/kg -389.5202
## 4: -393.5989 NA mgO2/hr/kg -393.5989
## 5: -390.0762 NA mgO2/hr/kg -390.0762
## 6: -365.5191 NA mgO2/hr/kg -365.5191
## -----------------------------------------Results
results <- data.frame(Clutch = Clutch,
Replicate =Replicate,
Male=Male,
Female=Female,
Population = Population,
Tank = Tank,
Mass = mass,
Chamber = chamber,
System = system1,
Volume = chamber_vol,
Date_tested = Date_tested,
Date_analysed =Date_analysed,
Swim = Swim,
Salinity = salinity,
Temperature = as.numeric(unique(firesting2$temperature)),
Resting_kg = mean(apoly_rmr$rate.output*-1),
Resting = mean(apoly_rmr$rate.output*-1)*mass,
rsqrest =mean(apoly_rmr$rsq))
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 59 | 6 | CARL354 | CARL355 | Arlington reef | 281 | 0.0006868 | ch3 | Asus | 0.04873 | 2023-04-04 | 2024-06-26 | good/good | 36 | 27 | 388.3666 | 0.2667302 | 0.9986 |
Maximum oxygen consumption
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch3
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.06 1.63
## -----------------------------------------
Subset data
cycle1.start <- Cycle_1.mmr[1,1]
cycle1.end <- tail(Cycle_1.mmr, n=1)[1,1]
cycle1.start.row <- which(firesting2_mmr$TIME == cycle1.start); cycle1.start## Warning in which(firesting2_mmr$TIME == cycle1.start): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="## Warning in which(firesting2_mmr$TIME == cycle1.end): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="cycle1_data <- firesting2_mmr |>
subset_data(from = cycle1.start.row,
to = cycle1.end.row,
by = "row") ## subset_data: Multi-column dataset detected in input!
## subset_data is generally intended to subset data already passed through inspect(), or 2-column data frames where time and oxygen are in columns 1 and 2 respectively.
## Subsetting will proceed anyway based on this assumption, but please ensure you understand what you are doing.## inspect: Applying column default of 'time = 1'## inspect: Applying column default of 'oxygen = 2'## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch3
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 2 4 5 6 7 10 11 12 13 14 15 16 17 19 20 26 27 28 32 33
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.06 1.20
## -----------------------------------------
Calculating MMR
## Warning: auto_rate: Multi-column dataset detected in input. Selecting first two columns by default.
## If these are not the intended data, inspect() or subset the data frame columns appropriately before running auto_rate()
##
## # summary.auto_rate # -------------------
##
## === Summary of Results by Highest Rate ===
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 421.6183 -0.07642660 0.9973620 NA 49 105 4238.44
## 2: NA 2 421.4139 -0.07637809 0.9973167 NA 50 106 4239.51
## 3: NA 3 421.2773 -0.07634764 0.9972568 NA 48 104 4237.37
## 4: NA 4 420.5607 -0.07617744 0.9970733 NA 51 107 4240.59
## 5: NA 5 419.8114 -0.07600108 0.9968315 NA 52 108 4241.66
## ---
## 220: NA 220 148.4036 -0.01366572 0.7675736 NA 220 276 4422.30
## 221: NA 221 144.9514 -0.01289049 0.7455492 NA 221 277 4423.37
## 222: NA 222 141.9209 -0.01221053 0.7321953 NA 222 278 4424.45
## 223: NA 223 139.4751 -0.01166202 0.7214833 NA 223 279 4425.53
## 224: NA 224 137.5559 -0.01123192 0.7139445 NA 224 280 4426.60
## endtime oxy endoxy rate
## 1: 4298.44 97.573 93.191 -0.07642660
## 2: 4299.51 97.565 93.114 -0.07637809
## 3: 4297.37 97.535 93.226 -0.07634764
## 4: 4300.59 97.501 93.105 -0.07617744
## 5: 4301.66 97.398 93.039 -0.07600108
## ---
## 220: 4482.30 88.171 87.388 -0.01366572
## 221: 4483.37 88.161 87.368 -0.01289049
## 222: 4484.45 88.149 87.297 -0.01221053
## 223: 4485.53 88.086 87.256 -0.01166202
## 224: 4486.60 88.036 87.210 -0.01123192
##
## Regressions : 224 | Results : 224 | Method : highest | Roll width : 60 | Roll type : time
## -----------------------------------------Adjusting
## adjust_rate: Rate adjustments applied using "mean" method.##
## # print.adjust_rate # -------------------
## NOTE: Consider the sign of the adjustment value when adjusting the rate.
##
## Adjustment was applied using the 'mean' method.
##
## Rank 1 of 224 adjusted rate(s):
## Rate : -0.0764266
## Adjustment : -0.0008082058
## Adjusted Rate : -0.0756184
##
## To see other results use 'pos' input.
## To see full results use summary().
## -----------------------------------------Converting units
mmr_adj2 <- mmr_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253)## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.selecting rates
mmr_final <- mmr_adj2 |>
select_rate(method = "rsq", n=c(0.93,1)) |>
select_rate(method = "highest", n=1) |>
plot(type="full") |>
summary(export=TRUE)## select_rate: Selecting rates with rsq values between 0.93 and 1...## ----- Selection complete. 14 rate(s) removed, 210 rate(s) remaining -----## select_rate: Selecting highest 1 *absolute* rate values...## ----- Selection complete. 209 rate(s) removed, 1 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time endtime
## 1: NA 1 421.6183 -0.0764266 0.997362 NA 49 105 4238.44 4298.44
## oxy endoxy rate adjustment rate.adjusted rate.input oxy.unit
## 1: 97.573 93.191 -0.0764266 -0.0008082058 -0.0756184 -0.0756184 %Air
## time.unit volume mass area S t P rate.abs rate.m.spec
## 1: sec 0.04873 0.0006868 NA 36 27 1.013253 -0.8611689 -1253.886
## rate.a.spec output.unit rate.output
## 1: NA mgO2/hr/kg -1253.886
## -----------------------------------------Results
results <- results |>
mutate(Max_kg = mmr_final$rate.output*-1,
Max = (mmr_final$rate.output*-1)*mass,
rsqmax =mmr_final$rsq,
AAS_kg = Max_kg - Resting_kg,
AAS = Max - Resting,
Notes=Notes,
True_resting="")
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest | Max_kg | Max | rsqmax | AAS_kg | AAS | Notes | True_resting |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 59 | 6 | CARL354 | CARL355 | Arlington reef | 281 | 0.0006868 | ch3 | Asus | 0.04873 | 2023-04-04 | 2024-06-26 | good/good | 36 | 27 | 388.3666 | 0.2667302 | 0.9986 | 1253.886 | 0.8611689 | 0.997362 | 865.5194 | 0.5944387 | ||
| ### Expor | ting data |
## Rows: 270 Columns: 25
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (9): Male, Female, Population, Chamber, System, Date_tested, Swim, Note...
## dbl (16): Clutch, Replicate, Tank, Mass, Volume, Date_analysed, Salinity, Te...
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.7
Enter specimen data
Replicate = 7
mass = 0.0007026
chamber = "ch2"
Swim = "good/good"
chamber_vol = chamber2_asus
system1 = "Asus"
Notes=""
##--- time of trail ---##
experiment_mmr_date_asus <- "04 April 2023 02 10PM/Oxygen"
experiment_mmr_date2_asus <- "04 April 2023 02 10PM/All"
firesting_mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_mmr_date_asus,"data raw/Firesting.txt"),
delim = "\t", escape_double = FALSE,
col_types = cols(`Time (HH:MM:SS)` = col_time(format = "%H:%M:%S"),
`Time (s)` = col_number(), Ch1...5 = col_number(),
Ch2...6 = col_number(), Ch3...7 = col_number(),
Ch4...8 = col_number()), trim_ws = TRUE,
skip = 19) ## New names:
## • `Ch1` -> `Ch1...5`
## • `Ch2` -> `Ch2...6`
## • `Ch3` -> `Ch3...7`
## • `Ch4` -> `Ch4...8`
## • `Ch 1` -> `Ch 1...9`
## • `Ch 2` -> `Ch 2...10`
## • `Ch 3` -> `Ch 3...11`
## • `Ch 4` -> `Ch 4...12`
## • `('C)` -> `('C)...15`
## • `('C)` -> `('C)...16`
## • `Ch 1` -> `Ch 1...18`
## • `Ch 2` -> `Ch 2...19`
## • `Ch 3` -> `Ch 3...20`
## • `Ch 4` -> `Ch 4...21`
## • `Ch1` -> `Ch1...22`
## • `Ch2` -> `Ch2...23`
## • `Ch3` -> `Ch3...24`
## • `Ch4` -> `Ch4...25`
## • `Ch1` -> `Ch1...26`
## • `Ch2` -> `Ch2...27`
## • `Ch3` -> `Ch3...28`
## • `Ch4` -> `Ch4...29`
## • `` -> `...31`## Warning: One or more parsing issues, call `problems()` on your data frame for details,
## e.g.:
## dat <- vroom(...)
## problems(dat)Cycle_1.mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_mmr_date2_asus,"slopes/Cycle_1.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE) ## New names:
## • `` -> `...8`Background rates
Pre-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",preexperiment_date_asus,"slopes"))
pre_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_pre1 <- pre_cycle1 %>% calc_rate.bg()##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] 3.358073e-05post-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",postexperiment_date_asus,"slopes"))
post_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_post1 <- post_cycle1 %>% calc_rate.bg() ##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.0005928561Resting metabolic rate
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch2
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.06 4.06
## -----------------------------------------
#### subset data
Tstart.row=which(firesting2_asus$TIME ==Cycle_1_asus$Time[1], firesting_asus$TIME)
Tstart.dTIME=as.numeric(firesting2_asus[Tstart.row, "dTIME"])
Tend.row=which(firesting2_asus$TIME ==tail(Cycle_last_asus$Time, n=1), firesting_asus$TIME)
Tend.dTIME=as.numeric(firesting2_asus[Tend.row, "dTIME"])
apoly_insp <- firesting2_asus |>
subset_data(from=Tstart.dTIME,
to=Tend.dTIME,
by="time")
apoly_insp <- inspect(apoly_insp, time=1, oxygen=2)##
## # print.inspect # -----------------------
## dTIME ch2
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 3 6 7 8 9 10 11 13 14 16 17 19 20 22 24 25 26 28 29
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.06 3.52
## -----------------------------------------
Extract rates
apoly_cr.int <- calc_rate.int(apoly_insp,
starts=(195+45+300),
wait=45,
measure=245,
by="time",
plot=TRUE) ##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from all replicates ...
## -----------------------------------------adjust rates for background
apoly_cr.int_adj <- adjust_rate(apoly_cr.int,
by = bg_pre,
by2 = bg_post,
time_by = Tstart.row,
time_by2 = Tend.row,
method = "linear")## Warning: adjust_rate: One or more of the timestamps for the rate(s) in 'x' do not lie between the timestamps for the 'by' and 'by2' background rates.
## Ensure this is correct. The adjustment value has been calculated regardless by extrapolating outside the background rates time window.## Warning: adjust_rate: background rates in 'by' and 'by2' differ in sign (i.e. one is +ve, one is -ve).
## Ensure this is correct. The 'linear' adjustment has been performed regardless.## adjust_rate: Rate adjustments applied using "linear" method.Converting units
apoly_cr.int_adj2 <- apoly_cr.int_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253) ## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.Plot curve
ggplot(as.data.frame(apoly_cr.int_adj2$summary), aes(x=row, y=rate.output*-1)) +
geom_point() +
stat_smooth(method = "lm", formula = y~poly(x, 2), color="red") +
theme_classic()
Rate filtering
apoly_rmr <- apoly_cr.int_adj2 |>
select_rate(method ="rsq", n=c(0.95,1)) |>
select_rate(method="lowest", n=6) |>
plot(type="full") |>
summary(export = TRUE)## select_rate: Selecting rates with rsq values between 0.95 and 1...## ----- Selection complete. 0 rate(s) removed, 20 rate(s) remaining -----## select_rate: Selecting lowest 6 *absolute* rate values...## ----- Selection complete. 14 rate(s) removed, 6 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: 10 1 349.1263 -0.01983227 0.962 NA 4569 4798 12609.66
## 2: 16 1 395.5926 -0.01870671 0.979 NA 7590 7819 15848.89
## 3: 17 1 430.5915 -0.02022344 0.993 NA 8094 8322 16389.30
## 4: 18 1 429.5327 -0.01951384 0.994 NA 8595 8824 16929.09
## 5: 19 1 431.6052 -0.01903446 0.979 NA 9099 9327 17469.40
## 6: 20 1 452.0224 -0.01958824 0.990 NA 9602 9831 18008.80
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 12855.08 99.258 94.855 -0.01983227 -0.0003208868 -0.01951138 -0.01951138
## 2: 16094.30 99.200 94.275 -0.01870671 -0.0005271242 -0.01817959 -0.01817959
## 3: 16633.72 99.130 94.318 -0.02022344 -0.0005614999 -0.01966194 -0.01966194
## 4: 17174.54 99.245 94.805 -0.01951384 -0.0005959005 -0.01891794 -0.01891794
## 5: 17713.85 99.137 94.707 -0.01903446 -0.0006302695 -0.01840419 -0.01840419
## 6: 18254.29 99.493 94.752 -0.01958824 -0.0006646455 -0.01892359 -0.01892359
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04958 0.0007026 NA 36 27 1.013253 -0.2260784
## 2: %Air sec 0.04958 0.0007026 NA 36 27 1.013253 -0.2106469
## 3: %Air sec 0.04958 0.0007026 NA 36 27 1.013253 -0.2278229
## 4: %Air sec 0.04958 0.0007026 NA 36 27 1.013253 -0.2192022
## 5: %Air sec 0.04958 0.0007026 NA 36 27 1.013253 -0.2132494
## 6: %Air sec 0.04958 0.0007026 NA 36 27 1.013253 -0.2192677
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -321.7740 NA mgO2/hr/kg -321.7740
## 2: -299.8105 NA mgO2/hr/kg -299.8105
## 3: -324.2569 NA mgO2/hr/kg -324.2569
## 4: -311.9872 NA mgO2/hr/kg -311.9872
## 5: -303.5146 NA mgO2/hr/kg -303.5146
## 6: -312.0804 NA mgO2/hr/kg -312.0804
## -----------------------------------------Results
results <- data.frame(Clutch = Clutch,
Replicate =Replicate,
Male=Male,
Female=Female,
Population = Population,
Tank = Tank,
Mass = mass,
Chamber = chamber,
System = system1,
Volume = chamber_vol,
Date_tested = Date_tested,
Date_analysed =Date_analysed,
Swim = Swim,
Salinity = salinity,
Temperature = as.numeric(unique(firesting2$temperature)),
Resting_kg = mean(apoly_rmr$rate.output*-1),
Resting = mean(apoly_rmr$rate.output*-1)*mass,
rsqrest =mean(apoly_rmr$rsq))
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 59 | 7 | CARL354 | CARL355 | Arlington reef | 281 | 0.0007026 | ch2 | Asus | 0.04958 | 2023-04-04 | 2024-06-26 | good/good | 36 | 27 | 314.7226 | 0.2211241 | 0.9836 |
Maximum oxygen consumption
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch2
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.06 4.06
## -----------------------------------------
Subset data
cycle1.start <- Cycle_1.mmr[1,1]
cycle1.end <- tail(Cycle_1.mmr, n=1)[1,1]
cycle1.start.row <- which(firesting2_mmr$TIME == cycle1.start); cycle1.start## Warning in which(firesting2_mmr$TIME == cycle1.start): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="## Warning in which(firesting2_mmr$TIME == cycle1.end): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="cycle1_data <- firesting2_mmr |>
subset_data(from = cycle1.start.row,
to = cycle1.end.row,
by = "row") ## subset_data: Multi-column dataset detected in input!
## subset_data is generally intended to subset data already passed through inspect(), or 2-column data frames where time and oxygen are in columns 1 and 2 respectively.
## Subsetting will proceed anyway based on this assumption, but please ensure you understand what you are doing.## inspect: Applying column default of 'time = 1'## inspect: Applying column default of 'oxygen = 2'## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch2
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 3 6 7 8 9 10 11 13 14 16 17 19 20 22 24 25 26 28 29
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.06 1.12
## -----------------------------------------
Calculating MMR
## Warning: auto_rate: Multi-column dataset detected in input. Selecting first two columns by default.
## If these are not the intended data, inspect() or subset the data frame columns appropriately before running auto_rate()
##
## # summary.auto_rate # -------------------
##
## === Summary of Results by Highest Rate ===
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 525.8170 -0.05509974 0.9930446 NA 58 114 7765.50
## 2: NA 2 525.8116 -0.05509827 0.9930375 NA 59 115 7766.57
## 3: NA 3 525.3794 -0.05504210 0.9928230 NA 60 116 7767.65
## 4: NA 4 524.3595 -0.05491049 0.9922750 NA 61 117 7768.72
## 5: NA 5 523.3051 -0.05477817 0.9919598 NA 57 113 7764.42
## ---
## 220: NA 220 295.8119 -0.02579178 0.9594312 NA 139 195 7852.55
## 221: NA 221 295.6749 -0.02577449 0.9594202 NA 140 196 7853.63
## 222: NA 222 295.2174 -0.02571764 0.9595988 NA 143 199 7856.86
## 223: NA 223 294.5886 -0.02563706 0.9599813 NA 141 197 7854.70
## 224: NA 224 293.8765 -0.02554721 0.9605398 NA 142 198 7855.78
## endtime oxy endoxy rate
## 1: 7825.50 97.710 94.689 -0.05509974
## 2: 7826.57 97.727 94.747 -0.05509827
## 3: 7827.65 97.691 94.706 -0.05504210
## 4: 7828.72 97.642 94.692 -0.05491049
## 5: 7824.42 97.694 94.749 -0.05477817
## ---
## 220: 7912.55 93.271 91.746 -0.02579178
## 221: 7913.63 93.282 91.689 -0.02577449
## 222: 7916.86 93.205 91.479 -0.02571764
## 223: 7914.70 93.339 91.641 -0.02563706
## 224: 7915.78 93.313 91.579 -0.02554721
##
## Regressions : 224 | Results : 224 | Method : highest | Roll width : 60 | Roll type : time
## -----------------------------------------Adjusting
## adjust_rate: Rate adjustments applied using "mean" method.##
## # print.adjust_rate # -------------------
## NOTE: Consider the sign of the adjustment value when adjusting the rate.
##
## Adjustment was applied using the 'mean' method.
##
## Rank 1 of 224 adjusted rate(s):
## Rate : -0.05509974
## Adjustment : 3.358073e-05
## Adjusted Rate : -0.05513332
##
## To see other results use 'pos' input.
## To see full results use summary().
## -----------------------------------------Converting units
mmr_adj2 <- mmr_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253)## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.selecting rates
mmr_final <- mmr_adj2 |>
select_rate(method = "rsq", n=c(0.93,1)) |>
select_rate(method = "highest", n=1) |>
plot(type="full") |>
summary(export=TRUE)## select_rate: Selecting rates with rsq values between 0.93 and 1...## ----- Selection complete. 12 rate(s) removed, 212 rate(s) remaining -----## select_rate: Selecting highest 1 *absolute* rate values...## ----- Selection complete. 211 rate(s) removed, 1 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 525.817 -0.05509974 0.9930446 NA 58 114 7765.5
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 7825.5 97.71 94.689 -0.05509974 3.358073e-05 -0.05513332 -0.05513332
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04958 0.0007026 NA 36 27 1.013253 -0.6388298
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -909.2369 NA mgO2/hr/kg -909.2369
## -----------------------------------------Results
results <- results |>
mutate(Max_kg = mmr_final$rate.output*-1,
Max = (mmr_final$rate.output*-1)*mass,
rsqmax =mmr_final$rsq,
AAS_kg = Max_kg - Resting_kg,
AAS = Max - Resting,
Notes=Notes,
True_resting="")
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest | Max_kg | Max | rsqmax | AAS_kg | AAS | Notes | True_resting |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 59 | 7 | CARL354 | CARL355 | Arlington reef | 281 | 0.0007026 | ch2 | Asus | 0.04958 | 2023-04-04 | 2024-06-26 | good/good | 36 | 27 | 314.7226 | 0.2211241 | 0.9836 | 909.2369 | 0.6388298 | 0.9930446 | 594.5142 | 0.4177057 | ||
| ### Expor | ting data |
## Rows: 271 Columns: 25
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (9): Male, Female, Population, Chamber, System, Date_tested, Swim, Note...
## dbl (16): Clutch, Replicate, Tank, Mass, Volume, Date_analysed, Salinity, Te...
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.8
Enter specimen data
Replicate = 8
mass = 0.0006180
chamber = "ch1"
Swim = "good/good"
chamber_vol = chamber1_asus
system1 = "Asus"
Notes="check max"
##--- time of trail ---##
experiment_mmr_date_asus <- "04 April 2023 01 42PM/Oxygen"
experiment_mmr_date2_asus <- "04 April 2023 01 42PM/All"
firesting_mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_mmr_date_asus,"data raw/Firesting.txt"),
delim = "\t", escape_double = FALSE,
col_types = cols(`Time (HH:MM:SS)` = col_time(format = "%H:%M:%S"),
`Time (s)` = col_number(), Ch1...5 = col_number(),
Ch2...6 = col_number(), Ch3...7 = col_number(),
Ch4...8 = col_number()), trim_ws = TRUE,
skip = 19) ## New names:
## • `Ch1` -> `Ch1...5`
## • `Ch2` -> `Ch2...6`
## • `Ch3` -> `Ch3...7`
## • `Ch4` -> `Ch4...8`
## • `Ch 1` -> `Ch 1...9`
## • `Ch 2` -> `Ch 2...10`
## • `Ch 3` -> `Ch 3...11`
## • `Ch 4` -> `Ch 4...12`
## • `('C)` -> `('C)...15`
## • `('C)` -> `('C)...16`
## • `Ch 1` -> `Ch 1...18`
## • `Ch 2` -> `Ch 2...19`
## • `Ch 3` -> `Ch 3...20`
## • `Ch 4` -> `Ch 4...21`
## • `Ch1` -> `Ch1...22`
## • `Ch2` -> `Ch2...23`
## • `Ch3` -> `Ch3...24`
## • `Ch4` -> `Ch4...25`
## • `Ch1` -> `Ch1...26`
## • `Ch2` -> `Ch2...27`
## • `Ch3` -> `Ch3...28`
## • `Ch4` -> `Ch4...29`
## • `` -> `...31`## Warning: One or more parsing issues, call `problems()` on your data frame for details,
## e.g.:
## dat <- vroom(...)
## problems(dat)Cycle_1.mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_mmr_date2_asus,"slopes/Cycle_1.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE) ## New names:
## • `` -> `...8`Background rates
Pre-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",preexperiment_date_asus,"slopes"))
pre_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_pre1 <- pre_cycle1 %>% calc_rate.bg()##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.0003237025post-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",postexperiment_date_asus,"slopes"))
post_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_post1 <- post_cycle1 %>% calc_rate.bg() ##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.0001686193Resting metabolic rate
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch1
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.06 4.06
## -----------------------------------------
#### subset data
Tstart.row=which(firesting2_asus$TIME ==Cycle_1_asus$Time[1], firesting_asus$TIME)
Tstart.dTIME=as.numeric(firesting2_asus[Tstart.row, "dTIME"])
Tend.row=which(firesting2_asus$TIME ==tail(Cycle_last_asus$Time, n=1), firesting_asus$TIME)
Tend.dTIME=as.numeric(firesting2_asus[Tend.row, "dTIME"])
apoly_insp <- firesting2_asus |>
subset_data(from=Tstart.dTIME,
to=Tend.dTIME,
by="time")
apoly_insp <- inspect(apoly_insp, time=1, oxygen=2)##
## # print.inspect # -----------------------
## dTIME ch1
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 3 6 7 8 9 10 11 13 14 16 17 19 20 22 24 25 26 28 29
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.06 3.52
## -----------------------------------------
Extract rates
apoly_cr.int <- calc_rate.int(apoly_insp,
starts=(195+45+300),
wait=45,
measure=245,
by="time",
plot=TRUE) ##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from all replicates ...
## -----------------------------------------adjust rates for background
apoly_cr.int_adj <- adjust_rate(apoly_cr.int,
by = bg_pre,
by2 = bg_post,
time_by = Tstart.row,
time_by2 = Tend.row,
method = "linear")## Warning: adjust_rate: One or more of the timestamps for the rate(s) in 'x' do not lie between the timestamps for the 'by' and 'by2' background rates.
## Ensure this is correct. The adjustment value has been calculated regardless by extrapolating outside the background rates time window.## adjust_rate: Rate adjustments applied using "linear" method.Converting units
apoly_cr.int_adj2 <- apoly_cr.int_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253) ## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.Plot curve
ggplot(as.data.frame(apoly_cr.int_adj2$summary), aes(x=row, y=rate.output*-1)) +
geom_point() +
stat_smooth(method = "lm", formula = y~poly(x, 2), color="red") +
theme_classic()
Rate filtering
apoly_rmr <- apoly_cr.int_adj2 |>
select_rate(method ="rsq", n=c(0.95,1)) |>
select_rate(method="lowest", n=6) |>
plot(type="full") |>
summary(export = TRUE)## select_rate: Selecting rates with rsq values between 0.95 and 1...## ----- Selection complete. 0 rate(s) removed, 20 rate(s) remaining -----## select_rate: Selecting lowest 6 *absolute* rate values...## ----- Selection complete. 14 rate(s) removed, 6 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: 5 1 189.2009 -0.009005217 0.986 NA 2054 2282 9908.80
## 2: 6 1 205.3678 -0.010069082 0.993 NA 2558 2786 10449.49
## 3: 8 1 226.0248 -0.010897361 0.986 NA 3563 3791 11529.23
## 4: 11 1 229.0694 -0.009792461 0.980 NA 5072 5301 13148.89
## 5: 15 1 258.2566 -0.010315282 0.981 NA 7087 7315 15309.11
## 6: 16 1 272.2180 -0.010835823 0.982 NA 7590 7819 15848.89
## endtime oxy endoxy rate adjustment rate.adjusted
## 1: 10153.34 99.845 97.698 -0.009005217 -0.0002785272 -0.008726690
## 2: 10694.03 99.972 97.662 -0.010069082 -0.0002700048 -0.009799077
## 3: 11773.69 100.190 97.518 -0.010897361 -0.0002529865 -0.010644374
## 4: 13394.28 100.100 97.676 -0.009792461 -0.0002274499 -0.009565011
## 5: 15553.84 100.090 97.597 -0.010315282 -0.0001934056 -0.010121877
## 6: 16094.30 100.200 97.657 -0.010835823 -0.0001848921 -0.010650931
## rate.input oxy.unit time.unit volume mass area S t P
## 1: -0.008726690 %Air sec 0.04565 0.000618 NA 36 27 1.013253
## 2: -0.009799077 %Air sec 0.04565 0.000618 NA 36 27 1.013253
## 3: -0.010644374 %Air sec 0.04565 0.000618 NA 36 27 1.013253
## 4: -0.009565011 %Air sec 0.04565 0.000618 NA 36 27 1.013253
## 5: -0.010121877 %Air sec 0.04565 0.000618 NA 36 27 1.013253
## 6: -0.010650931 %Air sec 0.04565 0.000618 NA 36 27 1.013253
## rate.abs rate.m.spec rate.a.spec output.unit rate.output
## 1: -0.0931011 -150.6490 NA mgO2/hr/kg -150.6490
## 2: -0.1045419 -169.1617 NA mgO2/hr/kg -169.1617
## 3: -0.1135600 -183.7540 NA mgO2/hr/kg -183.7540
## 4: -0.1020448 -165.1210 NA mgO2/hr/kg -165.1210
## 5: -0.1079857 -174.7342 NA mgO2/hr/kg -174.7342
## 6: -0.1136300 -183.8672 NA mgO2/hr/kg -183.8672
## -----------------------------------------Results
results <- data.frame(Clutch = Clutch,
Replicate =Replicate,
Male=Male,
Female=Female,
Population = Population,
Tank = Tank,
Mass = mass,
Chamber = chamber,
System = system1,
Volume = chamber_vol,
Date_tested = Date_tested,
Date_analysed =Date_analysed,
Swim = Swim,
Salinity = salinity,
Temperature = as.numeric(unique(firesting2$temperature)),
Resting_kg = mean(apoly_rmr$rate.output*-1),
Resting = mean(apoly_rmr$rate.output*-1)*mass,
rsqrest =mean(apoly_rmr$rsq))
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 59 | 8 | CARL354 | CARL355 | Arlington reef | 281 | 0.000618 | ch1 | Asus | 0.04565 | 2023-04-04 | 2024-06-26 | good/good | 36 | 27 | 175.3276 | 0.1083525 | 0.9844 |
Maximum oxygen consumption
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch1
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.06 1.63
## -----------------------------------------
Subset data
cycle1.start <- Cycle_1.mmr[1,1]
cycle1.end <- tail(Cycle_1.mmr, n=1)[1,1]
cycle1.start.row <- which(firesting2_mmr$TIME == cycle1.start); cycle1.start## Warning in which(firesting2_mmr$TIME == cycle1.start): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="## Warning in which(firesting2_mmr$TIME == cycle1.end): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="cycle1_data <- firesting2_mmr |>
subset_data(from = cycle1.start.row,
to = cycle1.end.row,
by = "row") ## subset_data: Multi-column dataset detected in input!
## subset_data is generally intended to subset data already passed through inspect(), or 2-column data frames where time and oxygen are in columns 1 and 2 respectively.
## Subsetting will proceed anyway based on this assumption, but please ensure you understand what you are doing.## inspect: Applying column default of 'time = 1'## inspect: Applying column default of 'oxygen = 2'## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch1
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 5 6 8 9 10 12 13 14 16 19 20 21 23 25 28 29 32 35
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.06 1.18
## -----------------------------------------
Calculating MMR
## Warning: auto_rate: Multi-column dataset detected in input. Selecting first two columns by default.
## If these are not the intended data, inspect() or subset the data frame columns appropriately before running auto_rate()
##
## # summary.auto_rate # -------------------
##
## === Summary of Results by Highest Rate ===
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 312.3218 -0.034494956 0.9976037 NA 224 280 6276.96
## 2: NA 2 311.5101 -0.034366130 0.9977515 NA 223 279 6275.89
## 3: NA 3 310.8862 -0.034267038 0.9978915 NA 222 278 6274.81
## 4: NA 4 310.3965 -0.034188509 0.9987325 NA 210 266 6261.96
## 5: NA 5 310.2106 -0.034158863 0.9987008 NA 211 267 6263.03
## ---
## 220: NA 220 147.7552 -0.007741457 0.8020670 NA 5 61 6041.90
## 221: NA 221 145.7775 -0.007415377 0.8035483 NA 4 60 6040.82
## 222: NA 222 143.3636 -0.007017374 0.8070180 NA 3 59 6039.75
## 223: NA 223 141.0573 -0.006636937 0.8186296 NA 2 58 6038.68
## 224: NA 224 139.2945 -0.006346014 0.8289034 NA 1 57 6037.60
## endtime oxy endoxy rate
## 1: 6336.96 95.784 93.672 -0.034494956
## 2: 6335.89 95.831 93.717 -0.034366130
## 3: 6334.81 95.875 93.735 -0.034267038
## 4: 6321.96 96.275 94.255 -0.034188509
## 5: 6323.03 96.265 94.248 -0.034158863
## ---
## 220: 6101.90 100.950 100.370 -0.007741457
## 221: 6100.82 100.930 100.370 -0.007415377
## 222: 6099.75 100.950 100.380 -0.007017374
## 223: 6098.68 100.970 100.430 -0.006636937
## 224: 6097.60 100.990 100.450 -0.006346014
##
## Regressions : 224 | Results : 224 | Method : highest | Roll width : 60 | Roll type : time
## -----------------------------------------Adjusting
## adjust_rate: Rate adjustments applied using "mean" method.##
## # print.adjust_rate # -------------------
## NOTE: Consider the sign of the adjustment value when adjusting the rate.
##
## Adjustment was applied using the 'mean' method.
##
## Rank 1 of 224 adjusted rate(s):
## Rate : -0.03449496
## Adjustment : -0.0003237025
## Adjusted Rate : -0.03417125
##
## To see other results use 'pos' input.
## To see full results use summary().
## -----------------------------------------Converting units
mmr_adj2 <- mmr_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253)## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.selecting rates
mmr_final <- mmr_adj2 |>
select_rate(method = "rsq", n=c(0.93,1)) |>
select_rate(method = "highest", n=1) |>
plot(type="full") |>
summary(export=TRUE)## select_rate: Selecting rates with rsq values between 0.93 and 1...## ----- Selection complete. 22 rate(s) removed, 202 rate(s) remaining -----## select_rate: Selecting highest 1 *absolute* rate values...## ----- Selection complete. 201 rate(s) removed, 1 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 312.3218 -0.03449496 0.9976037 NA 224 280 6276.96
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 6336.96 95.784 93.672 -0.03449496 -0.0003237025 -0.03417125 -0.03417125
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04565 0.000618 NA 36 27 1.013253 -0.3645576
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -589.899 NA mgO2/hr/kg -589.899
## -----------------------------------------Results
results <- results |>
mutate(Max_kg = mmr_final$rate.output*-1,
Max = (mmr_final$rate.output*-1)*mass,
rsqmax =mmr_final$rsq,
AAS_kg = Max_kg - Resting_kg,
AAS = Max - Resting,
Notes=Notes,
True_resting="")
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest | Max_kg | Max | rsqmax | AAS_kg | AAS | Notes | True_resting |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 59 | 8 | CARL354 | CARL355 | Arlington reef | 281 | 0.000618 | ch1 | Asus | 0.04565 | 2023-04-04 | 2024-06-26 | good/good | 36 | 27 | 175.3276 | 0.1083525 | 0.9844 | 589.899 | 0.3645576 | 0.9976037 | 414.5714 | 0.2562051 | check max | |
| ### Expor | ting data |
## Rows: 272 Columns: 25
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (9): Male, Female, Population, Chamber, System, Date_tested, Swim, Note...
## dbl (16): Clutch, Replicate, Tank, Mass, Volume, Date_analysed, Salinity, Te...
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.